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1. To find the potential at an internal point P of a thin ellipsoidal shell bounded 
by similar concentric ellipsoids, the law of force being the inverse /cth power of the 
distance. 

Let (x 9 y, z) be the coordinates of P ; (I, m, n) the direction cosines of a radius 
vector drawn from P to any point Q on the surface ; and let PQ = r. Let (a, 6, c) be 
the semi-axes of the ellipsoids, and {a, /3, y) their squared reciprocals. Then from the 
equation of the ellipsoid 

{aP + /3m 3 -f yn % ) r 3 + 2 (alx + j3my -j- ynz) r — ( 1 — ax 2 — $y % — yz 2 ) — . (1 ). 
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For the sake of brevity we shall put 



1/D 2 
F 
E 



,2 



aP + /3m 2 + yri* 
alx + fimy + ynz 



r • 



1 — ax % — /3y 2 — y2 2 



• (2). 



Thus D is the semi-diameter parallel to the radius vector r, F is a linear function of 
(x, y, z) such that F = represents a plane through the centre O diametral to 
PQ, and E = is the equation of the ellipsoid. 

The quadratic may therefore be written in either of the forms 



r 2 + 2FD 2 r - ED 2 = 



u 2 — 



2F 
E 



1 



u 



ED* 










* • 



(3) 



where u is the reciprocal of r. The roots r v r % have opposite signs and 



r x + r 



2 



- 2FD 2 , r x r 2 



ED ; 



(4). 



Taking together the two elements of the shell enclosed by an elementary cone 
whose vertex is the attracted point P and whose solid angle is do), the potential of 
the whole shell is evidently 



V = 



K 



1 

- 1. 



< 



p Y r^ day d?\ p 3 r 3 3 day ( — - dr 2 ) 



it K 



i-l'l)"- 1 



■ (5), 



where r l9 r 2 have their proper signs as given by the standard quadratic (3), r x being 
the positive root and r 3 the negative root. Also p l9 p 2 are the densities of the shell 
at the extremities of the chord through P. 

It is evident that if we integrate all round the point P every element of the 
surface is taken twice over. We must therefore halve the result given by (5). 

2. The values of dr v dr % depend on the kind of shell under consideration. 
Supposing that the shell is a thin homoeoid, we differentiate the quadratic (l) on the 
supposition that b/a and c/a are constant. We thus find 



dr = 



da 



a, 



I) 3 



' r + M3 3 ' 



Using the relations (2) and (4), it immediately follows that 



do\ 



da 2D 3 



dr: 



a q\ 



2 



aa 

a r, 



2D 2 



(6). 



2 



These values of dr v dr % may easily be verified by a short geometrical proof. 
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The potential at an internal point, P, of a thin homoeoid is therefore 



V 



1 da- 



te — 1 a 



r /Vi 3 ~ K + p*A-n?~ K 



r l ~ r 2 



D^dco. . '. . . . (7). 



where the integration extends all round the point P, Writing da/a = /x, this formula 
also represents the potential of a thin stratum placed on the ellipsoid whose surface- 
density at any point Q m'fipp, where p is the perpendicular from the centre on the 
tangent plane at Q. 

3. If the laio of force is any even poiver of the distance, k is even, and if the density 
is uniform and equal to unity the expression (7) takes the form 



/e — 1 a 



Z — K /y> 3 — K 

L. 13 

i\ - r a 



T)*da) (6), 



in which the subject of integration is a symmetrical function of the roots of the 
fundamental quadratic (3). Since such a function can be expressed as an integral 
rational function of the roots, the integration is much facilitated. 

If k be greater than 3, it will be more convenient to write r = 1/u. Supposing k 
to be an even integer the formula (7) then becomes 



_ _J_ da 1 

k — 1 a E 



Pl^l — Pofi'2 7 / rv \ 



u i ~~~ % h 



If the law of force is any odd power of the distance and the density is uniform, 
the potential at an internal point becomes 






da 



/> 



k — 1 a 



Z — K I nt Z — K 



n* 6 K _l n 



/}■» — Of 

'1 f % 



D*dco . , . . . , . (10). 



This is not a symmetrical function of the roots of the quadratic. It appears, 
therefore, that the odd powers of the force lead to a more complicated analysis than 
the even ones. 

4. We may apply a similar method to solid ellipsoids. Not to complicate matters 
at the beginning of the investigation, let us suppose that the ellipsoid is homogeneous 
and of unit density. The potential of the elementary cone at its vertex is found by 
integrating from r = to r. The potential of the whole solid at an internal point P 
is therefore 

V = i ( 73T^- 4 ) ff «-* + (~ *)<-) da, + C . . . . (11), 

where the result has been halved for the same reason as before. Writing 
do) = sin $ d6 d<j> the integration extends from <f> = to 2n 9 6 == to tt. 
mdcccxcv. — a. 5 z 
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When k is negative or positive and less than 4, the constant C is zero. When k 
is positive and greater than 4 the constant is infinite. Since the potential of a small 
sphere of fixed radius e and unit density at its centre is an absolute constant, we 
may take the lower limit in the integration for r to be e, and omit the potential of 
the small sphere. The value of C will then be finite and depend on the magnitude 
of e. As it will disappear when V is differentiated to find the forces, we shall omit 
it here also. 

5. Before discussing the general values of the expressions (7) and (8), it will be 
useful to notice a few particular cases. 

Suppose that the ellipsoidal shell is homogeneous and that the force varies as the 
inverse fourth power of the distance. We then have k = 4, and the potential at any 
internal point P is 

"XT yC6 ff— D 3 , 4tT/X 1 

3 J J r x r 2 3 E 

by using the standard relations (4). Here, as before, 

E = 1 — ax 2 — fiy* — yz % . 

Thus the level surfaces inside the homoeoid are quadrics similar and similarly 
situated to the bounding surface. This agrees with the result given by Townsend, 
in the 'Quarterly Journal/ vol. 12, 1873, "on the attraction of the ellipsoid for the 
law of inverse fourth power of the distance/' 

The potential V at an external point (x'y'z') follows immediately from that at an 
internal point. We have 

4?r/i cibc 1 



V' = 



ab'c' E' 



where E' may be written on either of the forms 
and 

a 2 + e 2 "T £3 + € S "T # + e 8 

Here a', b\ c' are the serni-axes of the confocal through the attracted point. 

The potential of a solid ellipsoid at an internal point may also be easily found for 
the same law of force. The integration for r, in equation (11), takes a logarithmic 
form when k = 4, and we have 

V = i £ \ da> (log r, + log (- »- 9 )} = i J do log (ED 3 ) 
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by using the standard forms (4). Integrating and omitting the constant, viz. : 
| dco log D s ; we have 

V = — log E. 

From this we easily find, by differentiating with regard to x, that the force at cm 
internal point resolved parallel to x is 

4c7T 1 x 



X = 



3 E a 3 



By Ivory's theorem it immediately follows that the force at an external point 
(x\ y\ z) is 

•yt _ 4tt ale x' 1 lx' 2 y n z^\" 1 

~ ~ T iw 'a^V\a 7l + v* + y± 

These values of X, X 7 agree with those given by Townsend in the memoir already 
referred to. 

We may deduce from the values of X 7 , Y', 71 the potential at any external point. 
Since 



we have 



Hence 







a 1 - 


x' 2 y' 9, z* 

\ + € * + J3 + 6 2 + c 3 + 6 3 — 


1, 




2x 


z=z « 


^'2 «/2 

1 * -1 


z't ' 


1 de* 


a? + e 3 


[ (a 3 + e 3 ) 3 ' (& 3 + e 3 ) 3 ' (c 3 


+ 6 3 ) 3 J 


dx f 






v , 27r abc 1 cfo 3 
— "" 3 a'Vd e 3 cfc' ' 







It follows by an easy integration that the potential of a homogeneous solid ellipsoid 
at an external point for the law of the inverse fourth is 

^j f _ 2tt f e2 abc cfrK 

We notice that the potential is constant if e 3 is constant ; thus the external level 
surfaces are confocals. 

Let p f be the perpendicular from the centre on a tangent plane to a confocal, then 
2p 7 dp' = d\. The force normal to the confocal is therefore 

"F — — — — — — 4?r abc 'fL 



dX dp' 3 a'Vd e* 

where a, &', c r are the semi-axes of the confocal and e 2 = a 2 ' — a 3 . For points on the 

5 z 2 
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same eonfocal the normal force is therefore proportional to the perpendicular on the 
tangent plane. 

The law of the inverse fourth power derives some importance from a theorem of 
Jellett's. Jellett showed that, if V K be the potential when the law of force is the 
inverse K th power of the distance, 



7-> '/•! \ 



7 



. 1 / d? d* d 

« + a — (jc + 1) ( K _ 2) [da? ^ dif T dz?J Kt 

By using this theorem, we may deduce the potential for any inverse integral power 
greater than the square, when the potential for the two laws of the inverse distance 
and the inverse fourth power are completely known. This theorem was communi- 
cated to the British Association at its meeting in Dublin in 1857, and briefly alluded 
to, rather than explained, in the report at page 3. But Townsend, in the ' Quarterly 
Journal,' vol. 16, 1879, has called attention to the theorem and used it to find the 
potentials of circular rings and plates, &c, for various laws of force. 

In this way the potential, when the force varies as the inverse sixth power, maybe 
derived from the potential for the inverse fourth, but as some of these potentials have 
the function we have called E in the denominator, while others have the still more 
complicated function a'&'c'E' in the same position, the labour will be found to become 
very great as we proceed to higher laws of force. 

We notice that the potentials of a homogeneous thin homceoid, at both external 
and internal points, have been expressed as algebraic functions of the coordinates of 
the attracted point. It follows from Jellett's theorem, that these potentials for all 
even higher laws of force can also be expressed free from all signs of integration. A 
similar argument will obviously apply to several other cases. 

Townsend also observes that Jellett's theorem remains true when the number of 
variables is extended beyond three. Thus, if 



V, = 2 ™ , r* = (x, - a Y f + (x % ~~ c h f + . . . + (x n - a n f 



we find by simple differentiation that 

6, Suppose the ellipsoidal shell to be homogeneous and the force to vary as the 
inverse sixth power of the distance. We then have k = 6, and the potential at any 
internal point P is by equation (8) 

V — -£- f f (r & JL rr 4- r %) -?^~ • 

V — e VI • 7 l' 3 T" > 2 ) _ 3 « 3 

O J J ' 1 ' 2 
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Using the formulae (4) this becomes 



V = -f JT(4F 3 D* + ED 3 ) —^ 



5 



"a: 



1 



— I (ulx + fi m y + ynzf dot + -^ (ap + /8m a + 7 n% ) d m 



W 



E 



Writing da = sin 6 dd dcf>, as usual, and putting cos 2 6 for P, m 3 , r? in turn, these 
integrals follow at once ; the terms with the products Im, inn, nl being obviously- 
zero. We therefore have 



V 



fl 4-7T 



o 



4 (a% 3 + /3 3 ?/ 3 + 7% 2 ) 



L(i 



oar 






« + j8 + 7 



y3^2 — ryZ^Y ' (1 — tf# 3 — /3;2/ 3 — 7# 3 ) 3 ^ 



The potential at an external point x' 9 y', % is therefore 



V 



477/, «6c [p'\* 
15~ aW 



1 1 1 



,/2 



a" ' 6' 3 ' c' 2 ' e 3 



4 -4P*(S+£ + 



a 



6 /6 






where a', &', c' are the semi-axes of the confocal through the attracted point, 
€ s = a /2 — a 3 , and p' is the perpendicular on the tangent plane at the attracted point. 
The potential of a solid homogeneous ellipsoid at an internal point for the inverse 
sixth power is easily seen to be 



V = — 



m 



15 



u 2 n' 2 + /3 3 ?/ + 7 3 £ 3 . a + /3 4- 7 
Zi . ^,„ ~ "j — 



E 3 






The potential at an external point is given later. 

7. Another interesting case is when the index k of the law of force is zero or 
negative. 

Taking k = so that the force is independent of the distance, we see from 
equation (8) that the potential of a homogeneous thin homoeoid at an internal 
point is 



V = 



p 



K 



(r-f + r i r & + r %) D s do*. 



Substituting (exactly as in the last article) from the standard forms (4) we see at 
once that the potential takes the form 

V = Aa? s + By 2 + Cs* + D 



where A, B, C, D are constants which depend on the axes of the ellipsoid. 
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In the same way when k = — 2, or the force varies as the square of the distance, 
we find that V 'is a quartio function of the coordinates of the attracted point 
containing only even powers. 

8. To find the potential of a thin homogeneous hommoid at an internal point when 
the force varies as any even inverse power of the distance greater than the square. 

Putting u = 1/r and using the formulae (4) 3 the potential (8) becomes 



A* 



' nr s - W*" 8 



k - 1 E 



V = -^-'- ^ -^ ^-da>, 



%b x — U c 



where u i9 u% are the roots of the quadratic 

?r — _ w — ^^To = 0. 

E ED 3 

Writing the quadratic in the form u? .— p^ — q = 0, we have by a known theorem 

^« + V = p* + np> - » 2 + . . . + n ^y^E^ P* " v ^ + - • • . 

where L (/) = 1.2.3.../! This series stops at the first negative power of p 3 so 
that the number of terms is \n + 1 or \ (n + 1 ) according as n is even or odd. 
We find by simply differentiating that 

- J *- = — — («," + «// 

u i ~ u a n dp 

to 

the differentiation being performed on the supposition that u x u % = — q is constant. 
We thus find that 

%h „ %h ~P +(n~Z)p q+ . .. + L(/) . L( „_2/-l)^ S" + --- 

In order to find V, we require the value of the integral 

r /2FV-4-2// 1 V 



ff 4 2 / / 7 (r/2F\*-*-v/ i y 



* One case of this equation is given by Sebrit, who uses it to find the expansion of x n — \jm n in powers 
of x -f 1/aJ, when that of a? M + l/a? w is known. Probably the equation for this case was suggested by differ- 
entiating cos nO. See ' Cours d'Algebre Superieure/ fourteenth lesson. 
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We know that 



(\l + [im + vnf* dco = (X 2 + /x 2 + z/ 2 ) 5 - 



47T 



2s + 1* 



Hence in our case 



\(alx + /3m;y + ynzY* do* = (a 2 ^ 2 + jS^y 3 + yV) s - -, 



Operating with 



V 



1 d 3 1 <P 1 # 

a dx~ /3 (hf 7 ch" 



on both sides, we have 



(odx + /3r/iy + ynzY s "^ (od 2 + i8m 2 + yn 2 ) d 
Repeating this process, and writing 25 = k — 4, 



47T 



(2s + l)2s(2s-l) 



[ 



F*-*- 2f ( ^ J" da> 



IV, 4ttL(«-4 -2/) 



L O - 3) 



V/P, 



where 



P = (A 3 + £y + y Vp- 4> . 



v V JL. • 



The potential of the thin homceoid when the force varies as the inverse nth power of 
the distance is therefore 



V 



9, 



TTfl 



'2\'- 



(k — 1) (* — 3) V K, 



1 E 

i I _i_ 

1 I 02 .. A I • 



E 2 V 3 



2 2 «-4 ' 2*1.2.(*-4)(* - 5) 



The general term within the bracket is 



1 L(*-4-/) 



2V L (/) . L (« - 4) 



WVf, 



"y" • » • 



P. 



and V, P have the meanings given above, viz., 



E = l - 



X/ 






y 



*2 






t) 



\ 






ce 



I 7.4. I „4, / 3 



# 



c 4 / 



v 



a dx^° dy* + ° d# 



It follows, therefore, that the potential of the homoeoid has been found free from 
all sign of integration. It has the form 



V=. A 



vp 



V 2 P 



1 



J*" 8 E"" E* ~ 5 * * * 
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where A, B, C are numerical factors whose values have just been found. The 
numerators of the powers of E are homogeneous functions of x 2 , y 2 , z* 9 whose 
dimensions are continually reduced by the differentiations implied in the operator V. 

The series has \ (k — 2) terms. Thus for the law of the inverse fourth it reduces 
to the one term given in Art. 5 ; for the law of the inverse sixth power it becomes 
the two terms found in Art. 6, and so on. 

The potential at an external point P\ whose coordinates are [x\ ?/, z'), follows from 
that at an internal point by what is practically a method of inversion. We describe 
through P' a confocal homoeoid whose axes are a\ b'\ c\ and of the same volume as 
the original homoeoid. The potentials of the tw r o homoeoids at corresponding points 
are then equal. Remembering the transformation already made in Art. 5, we find 
that the potential at an external point is 

V' — 2irfJb ab ° I -V" 3 Jl 4- — — 4-~ E /2 V /2 lp. 

~ fa _ 1)(k- 3) flW \E7 \ + 2 a *~4 + 2 4 1.2 (k -4)(* - 5) + * • 'J ' 

where E' and P' have either of the forms 

■^ — 1 "" a /4 &V " c '4 — ^ ^/4 + //4 + c '4 ) » 



V = 



« /6 ' U & l c' 6 

O 7/0 I 70 7/3 1^ 7 /O 



It should be noticed that the differentiations implied in the operator V' are to be 
performed on (x, y\ z f ) on the supposition that a\ b\ c are constant. 

9. To find the potential of a thin homogeneous homoeoid at an internal point when 
the force varies as the inverse Kth power of the distance and k = 2, i^ero, 07* any even 
negative number. 

We now use the potential (8) in the form 



V 






where r l9 r<> are the roots of the quadratic 

r* + 2FDV - ED 3 = 0. 
Puttings =■— 2FD 3 , g = ED 2 , we substitute in the series 

r ®~~ K t 3-/c ( '■ //if ■ //* 1 ^ 

-V-t- = ^ 2_K + (1 - *) p ' Kq + L if- i2 p - """ « + • • • 
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This series has ^ (4 — #) terms, and the general term has been written at length in 
Art. 8. 

To determine the potential we must therefore find the value of 

\\f~- K ~ 2f q f Wda> == W |j , (2F) 2 "«^(D 3 ) 3 — f dco 



where / has any positive value from zero to | {2 — k). To effect this integration we 
shall use the following lemma. 

10. Lemma. By a known theorem we have 



d 



CO 



oo _ t 



where 



Now 



aP -f /3m 3 -f <yn 2 



2lT 



*>o 



'v~*dv 

Q~ 5 



Q 3 = (a + v) (ft + v) (y + v). 



— J- C l~ JL ^ \ ~ ^ ^ 

cZa d£ dry/ Q ~~ <fe Q 



Hence, operating t times on the equation with the left-hand operator, we find 



d 



CO 



(ctP -f- /3m 3 + ryn?y +1 



2tt ( — 1)* 



,00 



/dV 1 
o • \^/ Q 



Differentiating again / times with regard to a, g times with regard to ft h times with 
regard to y, we have 



ffim*0n u day Tr f °° x _ / rf V 1 



where 






^TrJ-iy L (2/) L(2y) L(2ft) 1 
L (s + I) ' L (/) • L (</)" L (h) 2 2 * ' 



where s =/+ g + h and L (/) = 1 . 2 . 3 . . ./ 
By an expansion we deduce from this 



(XI + /xm + to) 38 dm 



(»P + 0m* + r / n i ) s+t + 1 



A 00 

= K' j v"~ h dv 



d\t 
dv, 



■HMtli! 



fjj 






V" 



_ a + # /3 + ^ <y + ?? 



* 1 



where 



K' 



Ifl2s) 2tt (- 1)' 



L (s) L (5 + 2 3 * 



The expansion is not difficult, for the odd powers of I, m, n, on the left hand side 

MDCCCXCV. — A. Q A 
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give zero on integration, and the genera] integral for any set of even powers is written 
down just above. 

The subject of integration on the right hand side is infinite when v = 0, but it is 
not difficult to see that the integral from zero to any very small value of v is 
ultimately zero. 

11. To apply the lemma to the integral required in our proposition we refer to 
Art. 1 for the definitions of F and D 2 . We have 

F = alx -f- /3my + ynz. 

We therefore put X = ax, p = fiy, v = yz, 2s = 2 — k — - 2f, s + t + 1 = 3 — k —f. 
This gives t = \ (2 — k) so that t is a positive integer when k is even and not greater 
than 2. We thus have 

V = ^i if * L (/) . V-V- -2/) • ^ W-^(D»)--^.. 

The summation extends from f= to \ (2 — #c) and .9 varies from -J (2 — k) to 0. 
The 3$ is therefore equivalent to an expansion by the binomial theorem. We thus 
find for the required potential 



* - 1 L (0 J W I ^« + ^/3 + *>Y + »J 

Substituting for E its value (given in Art. l), the potential becomes 












3 



where £ = -| (2 — k). 

Writing v — 1/u we find that the potential of a thin homoeoid when the index k of 
the force is even and not greater than 2 at an internal point, is 



d V abcu* T\f 



a: — 1 Jo V duj Q 1 L (0 

where 



o O O 

•y»,J o/" £-" / 

R = 1 — — 

a 2 + % 6 3 -f w c 3 + '# 

Qj* = (a 2 + w) (6* + w) (c 3 + u). 



At an external point (x'y'z') the potential is 
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V 



'7T/J, 



fC 






.CO 



(u — e 2 ) f du \ (u 



€*)* 



d 
du 



1 abc (u — € 2 )"^ / 



u 



f K'« 



Qi 



\u - e 2 / L (t) 



where B/ differs from R only in having (x f y f z) written for xyz. 

Here t = -| (2 — k), so that when £ = the formulae reduce to the well-known 
results for the law of the inverse square. 

It is easy to verify these results in the case of a spherical shell. 

12. To find the potential at an internal point P of a thin heterogeneous homceoid 
when the law of force is an inverse even positive power of the distance, and the density 
at any point (£ ??, £), is tf> (£ rj, £). 

The coordinates of the attracted point P being (as, y, z) 3 the coordinates of any point 
Q on the surface are f == x + lr 9 &c. If />j, p 2 are the densities at the extremities of 
a chord through P, we have 



Pi 






1.2 



where S = Z d/cfas + m d/cfa/ + « d/dz. In the same way 



P2 






JL . Zi 



The potential at P (as in Art. 3) is 






fi 1 



f /• « - 3 k -~ 3 



™1 - «8 



Substituting for ^ 1? /) 3 , this becomes 



V = 



M 



1 



# 



1 EJ 



?&■ 



^n 



* — 3 « - 3 k — 4 



u. 



It 



u. 






to, 



L (^ - 4) ^ 



+ zero + 



S*~ ^ 



L (k — 2) r x 






S* ~ ^ 






^•2 



i~ *W + LOT- 1) "*•!-" r 9 ^~ '™ + * • ' 



(i(W. 



The series is divided into two parts, one with w's and the other with r J s ? so that all 
the powers used are positive. The two parts will be represented by V 2 and V 2 , so 
that V = V x + V a . 

13. Ze£ us examine first the terms which contain u only. These may be written in 
the form 



^i= 



/* 



K 



1 E J 



i 



g/c-3-wJj B K "" 4: d> 

day \ <j&S K _ 3 + S^>S K _ 4 + . . . + L / ~g ™ x S» + • • • + l7^L"4) ^i 



where, as in Art. 8, 



6 A 2 
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S. = p-i + (n - 2) p*-*q + . . . + ^J±^ p - V-y + . . . 

P = ^ = ^7 (aa^ + #ym + yaw), <? = EJ)3 • 

Before attempting the general integration it is useful to examine what kind of 
integral we shall arrive at. We notice that every term of S K "" 3 "" tt ^.S» is an integral 
rational function of the direction cosines (I, m, n) of the chord. All the integra- 
tions can, therefore, be completely effected. The terms with odd powers of (I, m, n) 
give zero on integration, while those with even powers may be found by the theorem 

r (/ + <7 + a + I) 

Since every term of V x is of the same dimensions in (£, m, n), viz.. k — 4 dimen- 
sions, we notice that the denominator of this expression, viz., V (k — 4 + f ), is the 
same for every term. 

Supposing the density <j> (£ tj, Q "to be a homogeneous function of (g, 77, £) of 
i dimensions, we see that Y x will take the general form 



k — 4 4- i • ±l k — 6 + 'i i ■*■■*■ k — 8 + i 



i 1?k - 3 • "17* - 4 » 17 k - 5 « # * * 



where H, is a homogeneous function of (x, y, z) of s dimensions. If i > k — 4, 
the series ends when the exponent of E is unity, if i < k — 4 the series ends when 
the numerator reduces to a constant. 

Referring again to the general expression for the potential v, we see that where 
the highest power in <£ (as, y y z) is less than k — 2, the only terms in the potential are 
those with u. It follows that when the law of force is an even inverse K th power of 
the distance (k > 2), the potential of a heterogeneous thin homceoid can be found free 
from signs of integration, provided the law of density at (£, rj, £) is a function of 
(£, rj, Q of less than k — 2 dimensions. 

14. Let us now consider the general integration of the expression for V P The 
following lemma will be required. 

Beginning (as in Art. 8) with the equation 

(axl + /Sum + yznf 8 dco = (otV + /3^ + yV)* --- — •— 

w 1 Te the integration extends all round the origin, let 

1 <p 1 d? 1 cP 

a ^ 2 B dv % 7 f^ 2 

„, A d , B rf , cl 

V — ~r~ -+■ TT 7 -T" 7 

a dx /3 ^,V 7 dz 

where A, B, C are three arbitrary constants. 
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We then find 

j J h (,ZS — ^y — £/) 

_____ f f& (c*?x° + @y + 7 V) * 

— L (25 -f 1) 

where, as before L (s) = 1.2.3 ... 5. In this equation all the exponents are to be 

positive. 

d d d 
To apply this lemma we let A, B, C signify — , - , —-, which operate solely on the 

(mo Ct'U Cv% 

<j) and not on the other x, y, z which occur in the equation. The equation may then 
be written 

W-W-ff 1 1 \/ , „.„, (A 2 + £V + 7W 

I — \ §•? do> = 47r V^ V ^ ^_-_Z_jl__z. 

L(28-2/-^)VDV L (2s +1) 

We see by what precedes that it will be convenient to arrange the expression for 
Y l in a series of inverse powers of E. Let us, therefore, seek the term containing 

— ) . Putting n—f—l = h, we notice that =^~ — ^ r contributes the 

one term 



_ — 



77/. ? 



L(ra - A - 1) L(2A - n + 1) \D 2 / L (* - 3 - w) E 7 < 

where every exponent is positive. To find all the terms with the given power of E 
we sum this expression, beginning at .n = h + 1 and ending at the lesser of the two 
n = 2h + 1 or n = k — 3. 

Using the lemma, the corresponding term of V x becomes 

2tt^ , (jA)""^ 1 A /K ~ 3 ~ n (tt^ + flY + T 8 * 8 )**'" 4 * / AV + ' 

k-1 ^ ^L(w-A-l) L(*-3-7i) L(«-3) " \E/ 

where 2 implies a summation with regard to n. 

If 2-A + 1 > ac — 3 the summation extends from n = h + 1 to k — 3 and the 
several terms form a complete binomial series. If 2h + 1 < k — 3 the summation 
extends from n = h -\- I to 2h -{- 1. In this case the missing operators required to 
complete the binomial series are A 7if l _>_'* ~ 5 ~* 2 \ A 7M " 2 a'*"~ 6 ~ 2 \ &c. 

Since A, A' are differential operators of the orders d 2 /dx* and d/dx 9 while the 
subject is of the order x K ~ 4 , the missing operators reduce the subject to zero. We 
may, therefore, suppose the binomial series to be complete. 

It immediately follows that 
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where (as in Art. 8), 

and the £ implies summation from h = to h = k — 4. 

When the ellipsoidal shell is homogeneous we omit the operator A', and put the 
density <j£> equal to unity. The result then agrees with that already obtained in 
Art. 8. 

The differentiations with regard to x 9 y 9 z 9 which are implied by V, V', operate only 
on P, while the differentiations, A, B, C, in V operate only on <f>. When written at 
length the operator takes the form, 

(i V + V')- P * = * (i V). + . (i * I + 1 1 £ + . . .) ( J V,- P 

+ s — I- d -i — + + - 2 - -**._*- + \uvv-'P + 



a 3 ffo? 3 tfe 3 a/3 & dy dx dy , 



15. We may write the result in a different form so as to avoid having two kinds of 
operators. We have 

^^=i(i+^+Ki +B ) ,+ i(s+ o )'-T-l-~- 

Let V* represent the total operator V when the differentiations implied act on both 

P and <f>. Let V^ be the partial operator V when the differentiations act on <j> only, 

then 

V+ 2V'== V/ - V^. 

The potential of the heterogeneous homoeoid therefore becomes 

16. As an example, let us find the potential of a thin homoeoid at an internal point 
(x 9 y 9 z) when the force varies as the inverse fourth power and the density at (£ ??, £) 
is m£ 

The series has only one term given by h = '0, since k = 4. The result is seen at 
a glance to be 

V= — : r - T ,a. 



» E 



T/V&s the potential is the same as if the homoeoid were homogeneous and its density 
were mx. 

This curious result may be verified by independent reasoning. Let P be the given 
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internal point, PQ any radius vector of the shell. Let the direction cosines of PQ be 
(I, m, n ) and its length be r. The density at Q is therefore proportional to (x + Ir). 
The potential is 

v _ -jj-^ ■+— -JJ— r -• 

The first term is the potential of the homogeneous shell. The second is the x 
attraction of a homogeneous homoeoid where the law of force is the inverse square and 
is therefore zero. 

As another example, let us find the potential of a heterogeneous thin homoeoid 
when the force varies as the inverse sixth power and the density at (£, 77, £) is ml;*. 

The series has three terms given by h = 0, 1, 2. Hence, since k = 6, 

Now VP = 2 (a +/3 + y), V 2 P= 0, WP = 0, V'P = 2aA;r + &a, V' 3 P = 2A 3 + &c. 
Remembering that <j> = -m^ and A<jf> = dfyjdx, A 3 <£ = dP^/dx 2 , we find, after an 
easy reduction, 

~ "~15~~ [ E + ~~E3 "* E* J * 

It immediately follows that the potential at an external point x\ y\ z f for the same 
law of density is 

y_ 4tlT J x/mj j_i*L <A 1 A 1 /A jl. 1 4. A\ ^! L 1 A^l! 1 ^V 3 1 c% ' 3 \ ^ V2 A I 

~ " 15 " a?Vc' cfi J E' + \a'* &'*" c ' 3 / «' 3 E'* + \ a' fl + l*" + 7 6 "/ ^ E*} ' 



where 



T?' f%(— Jl. *. L - J — » 



and a', 6', c are the semi-axes of the confocal through the attracted point. 

17. As a third example, we may apply the formula to find the potential of a thin 
homogeneous focaloid when k is greater than 4. 

The thickness of a thin homoeoid is /xp, while that of a thin focaloid varies inversely 
as p. A homoeoid, whose density is ra/p 3 , is equivalent to a homogeneous focaloid 
of unit density whose ada = m/x. Putting then 

<f> = m/p* == m (aV + /3 2 ?/ 3 + y V), 

we have, if n = k — - 4 — h< 
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(V, - V + )» P<£ = { V/< - nVr 1 V,} P<£ 

= V" (P<£) - 2« (« + /3 + y) mV n ~ l P 

/ 1 \iU-2) I 1 U(«-4) 

= ™V (p) - 2n (a + + y ) mV- 1 (^ J 

Substituting this in the formula for V, already written down, we have the 
required potential. 

Suppose the law of force to be the inverse sixth power, the potential of a 
homogeneous thin focaloid, at an internal point, is 

v = ^ {.i±to! + 1 >{±mj + 4 (aW + w + /2i) ) + ±; ( 

and at an external point 

^, 4:7ram ctbo 1/1 , 1 , 1 , 4\ /0 

V = — - — • — -+-—-+- — 4~ — r n z 

15 aW e 4 \a* ^ 6 3 T c 3 e 3 / ^ ? 

where a', 6', c' are the axes of the confocal through the external point, and e s = a! % — a 2 . 
It appears that the potential varies as the square of the perpendicular p' on the 
tangent plane for all points on the same confocal 

18. When the density <j> (£ 77, £) of the thin homoeoid is a function of the 
coordinates of a degree higher than k — 2, the potential includes some of the terms 
of Art. 12, which contain r l and r g . We shall now examine these terms. In this 
way we shall arrive at an expression for the potential of a thin heterogeneous homoeoid 
at an internal point where the index k of the force is a positive even integer greater 
than the square. 

These additional terms are given by 



7 fJb I O <& o 1 O "(P 01 1 l O 



2 



A- — JL , 



<P ft 4- b "9 « rt 4- 4- __A__2 S + 

L (/c — 2) L L (« — 1) L (/c 4 71 — 3) 



^ ) c/g>. 



Ill 



S, = ^^ = p"-' + (n - 2)^-»g + . . . + r ,^7^uf^-¥+ • • • 



7\ - T % 



L(/).L(n-2/- 1) 

_p = — 2D 3 F, <7 = — r i r s ~ z ED2 ' 

and D, E, F are expressed as functions of the direction cosines (I, m, n) in Art. 1. 
Before attempting the general integration, it is convenient to examine the kind of 

integral we shall arrive at. 

Consider first the dimensions of the several terms of V a as regards the coordinates 
(as, y, z) of the attracted point. Let 4> (x, y, z) be homogeneous and of i dimensions, 
then remembering that F is a linear function of (x, y, z) we see that, if the terms are 
arranged in powers of q or E, V 2 takes the form 

V a = H,_ K+2 + KH W + E 2 H ; _ 2 + • • • 
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where H* is a homogeneous function of (as, y, z) of s dimensions. The series 
terminates when H, reduces to a constant. 

To calculate the values of H^ we must consider the dimensions of V 2 with regard 
to (Z, m, n). Since positive powers of D now occur in the series, the integrations 
cannot be completely effected. The coefficients of the powers of (as, y, z) in H, 
will be definite integrals of the same kind as those which occur when the force 
varies as the inverse square. 

The terms of V., which contain any odd exponents of (I, m, n) give zero after 
integration. Omitting these it will be seen that every term of V 2 is of the form 
pf m ty n M J)%(f+ </ + &)- (*~ty^ where jT, g, h here stand for any positive integers. Now 
putting t = \ (k — 2), the integral of any one of these terms may be found by the 
equation 



IT Wm frn** dm _ / d V/_d V / d \ 

]](al 2 + 13™*"+^ + *-' ~ \d*J \d/3j [dry) 



h f °° v ^-4 fa) 

o ~Q~~ 



where 



L (2t) L (* - *) ' J -ry-r r h 

Q 3 = (a + v) (/J + v) (y + v). 

In any particular example, therefore, the integrations may be effected (as far as 
they can be clone) by using this equation. 

1 9. To obtain the general integral of the series for V 2 , we require a lemma. 
We deduce from the integral given in the last article that 



i. 



where 



(*Z 3 + ^T^ 2 r 1 "' — Jo W + tf + /8 + tf + 7 + i>/ Q ' 

~ = L ( S ) L (s - 2*-* L (20 ' 
Q 2 = (a + v) (/3 + v) (y + v), 

and « is a positive integer less than s (see the authors 'Statics/ 1892, vol. 2, 
Arts. 219 and 220 for this integral and the one used in the last article), 
Let 

where A, B, C are any constants, we find 

([(^L^l P™ + vn )**~ 9 (A / + Bm + Cn)ff dco 
J J (*P + /3m 3 +^ 3 y+ T=r ^ " 



L(2 S -^)_. 4 rr x 8 , 

L (Js) Jo [a + v 



5 v i ~ i dv 
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To apply this lemma we suppose, as before, that A, B, C represent differentiations 
with regard to x, y, z, but to operate only on <£. Hence supposing X, fi, v to stand 
for ax, fiy, yz } and using the abbreviations given in Art. 1 5 we write the lemma in 
the form 



f. 






< + &C 

a + V J 



\ s v e " 



civ 



Q 



20. By substituting for S„ its value we find 



f 



>k + n — Z 



* 



L (# -f~ % -— 3) 



S* D 3 c^ou 



8 * + * '* ^/- r °! / ;- 1 V x -(-2F)»-^- l D' (n ^>E / . 



L (* + n - 3) ^ L/. L (w - 2/ - 1) 



To use the lemma we put # = k + n — 3, 2s — g = n — 2f 
These give 4 = 4 • (k — 2) so that k must not be less than 2. 
positive. The double integral therefore becomes 



They also make s — t 



(- 1)" 



-1 



27T L (0 H + ^ _ j 



L<> + 2*- 1) 2"- 1 L(2t) 






/ >3 



X' 



5 qjt ~~ -2 rim 



n 



f+t-1. 



where 2 (t + 1) has been written for k: and 5 = 

Since £ here implies summation from/= Q to f=^n— I or % (n — 1), the terms 
which follow the 2 do not form a complete binomial expansion. But, since the 
operator & n + 2t ~ l contains differential coefficients of the order (d/dk) n + 2t ~\ while the 
subject is of the order X 8(tt " /+ *"" 1) , the missing terms are zero. We may therefore 
sum the series from/ = to n + t — 1. We therefore have 






27T/JL L (I) 



-co 






u+2i-l 



R»+<-l 



«-~l L(20 J L(^ + 2tf-l) L(f6 + *- 1) 






where 



Z/" 



= Mi -r ~~~~ T /-> . i 

# -f w p -f- # 7 + ^ 



and 2 implies summation for integer values of w, beginning at n — 1 onwards until 
the terms vanish by the operation of the differential coefficients A, B 5 C, contained 
in A on the subject <f> (x, y, z). 

21. We may simplify this series by performing the operation a* + 2 *~ 1 R^"" 1 once 
for all. For the sake of brevity let 



S 



/ AX B^ Cv 

V a + v p + v y + Vj 



m 






-|- 



B 3 



a -\~ V ft -h V 






C 



^3 



7 



v 



We then find, when g < h, 



A'K* 



E ; *™^ 



+ 



E^ +1 S^~ 2 T 



+ 



I#-*+3S*MT 3 



L(y)."L(A) — L(A-^).L(^) nr L(A-flr + l).L0?-2) L(A- # + 2).L(#-4).L(2) 



+ &c, 
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the series stopping when the exponent of S becomes negative. The form of the 
expansion having been found, the simplest proof is that supplied by the method of 
induction. Performing A on both sides of this equation and noticing that 

AB^S* = mR^S^ 1 + nU m B n ^ 1 . 2T, 

we may see that since the theorem is true for g = 0, it is true for all positive integer 
values of g. It is unnecessary to reproduce the work. 

When g > h 9 some of the indices of R at the beginning of the series are negative. 
These terms are to be omitted and we then have 



&W 



g2A-^Qy-A 






Rg 2A-0-2'jy-A+l 



+ 



JpggA-^rjy-A+g 



L(g)L(h) L(2h-g).L(g-~h) ' L(2h-g-2)L(g-h + l) ' L(2)L(2A-if-4)L(//-A+2) 



+ &c. 



We thus have 



AB+2/-1 'Jtyl+t— 1 



L (n + 2t — 1) L (n + t — 1) 



_ 1 V" 1 _ (-jS)*" 1 ^ , R ( - |S)^ 3 F +1 1 



2 1 






L(n-l)L(t) l L(ra-3)L(* + 1) 2 8 



— |— ^ s 



In this equality S and T contain A, B, C, which represent differential coefficients 
of the function cf> (x, y y z). Also to find V 2 we have to sum both sides of this equality 
for all values of n from n = 1 onwards, until the terms become zero by the continued 
operation of S and T on <£, Since the series stops at the first negative power of S, 
all negative exponents are to be omitted. We therefore write the summation in the 
form 

A L(^ + 2^-l)L( ra + <-l) 1 *' ^T^)* L(»-l) t L(« + l)2» 2 ' L(»-3)~ "*"••• 

We notice that on the right-hand side the first t implies summation from n = 1 
to oo , the second from n = 3 to oo , and so on. Now by Taylor's theorem 






L (A) y \ a + V 



y 



v* 



P ~f v 



z 



v 



7 + i 



7 



where the 2 implies summation from h = to oo . We therefore have 



3 ""«-lL (it) J o ~"~Q~ 



T' 






RTHl 



ET< 



1 



L(t) ' L(* + 1) 2 3 



™j ... 



</>(x 



a -f V 



, &C ) . 



Here T is an operator on (#, y, 2). The quantities X, /jl, v are constants for which 
we write ax 9 fiy. yz after the operations indicated by T have been performed. We 
may reverse the order of these operations and write 

6 B 2 



918 DR. E. J, ROTTTH ON" THE ATTRACTION OP ELLIPSOIDS FOR 

X vx a vu 

/y> mmm , .... _. — — qt ^^ 1 . — „ v_ — Xrr* 

a + V a -f v ^ p -f V /3 + V 

before the operation T has been performed, provided we replace the d/dx, d/dy, d/dz 

-, t -, . m , a 4- V d ,8 -f v d <y + v d 

included in 1 by — , 7- , — respectively. 

J v ax v ay v dz l 

The operator T then takes the form 



IT, . v d* . , n . , d? . , . , (P 1 



T = - < (a + <y) To + (£ + t>) 7 1 ; + (7 + v) — [ r ; 

w l ^ /l " v ri vr ^ j 

also by the same substitution 

a 4- ?; /3 4- ^ 7 4- ^ 

which must now be placed on the left-hand side of the operator to indicate that it is 
not included in the subject. 

We thus find for the partial potential of a thin heterogeneous homoeoid, whose den- 
sity at the point (£ 17, £) is </> ('£ r), £), 

2tt/a_ L(0 fV-ta? f T^ RT^ + 1 !__ 1 / vx 

3 ~~ k ~ 1 L(20 ! " Q t l7(V) + L(* + 1) 2 3 + ' # ' J ^ W 4-T ? " ' 

To this we add the partial potential V l5 already found in Art, 14, for the complete 
potential. 

22. It is sometimes more convenient to use the axes a, h, c of the ellipsoid instead 
of the squared reciprocals a, /3, y. Let us also suppose that the density at any point 
(£ V> ' lB given by \jj (g/a, r)/b, Ijc), so that we may the more easily deduce the 
potential at an external point. Also writing v = 1 /u, we have 






Tp* ~« ^ . / 



a 3 + u W 1 4- ^ c? 4- ?6 ' 
T _ a?__+_u <P ¥ 4- v, d? <?_+_n d* 

Q^ = (cfi + tt) (& 3 +• w) (c 3 + it). 

Tfte partial potential of a thin heterogeneous homoeoid at an internal point is then 

V - 277>l L(£) pafo^ nyr , / g£ _&y__ _««_ 
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ivhere the operator M is 

— L (t) + L (t + 1) 2 s " + ' ' ' + L(n) . MtXn) '¥> + ' * * ' 

and! t = -|(/c — 2) ivhere k is the index of the law of force. 

Having found the law of the potential of a homoeoid for the inverse /c th power, 
the potential for the inverse (#c + 2 ) th power may be obtained by Jellett's theorem, 
already mentioned in Art. 5. We might thus obtain a new proof of the formula by 
the method of induction. The theorem has been verified in this way, but the 
demonstration is long and rather tedious. 

23. To deduce the potential of the heterogeneous homoeoid at an external point 
(as', y', z') 9 we follow the same process as before. Let a, b\ c be the axes of the 
confocal through the attracted point. Let a' 3 - a 3 = e 3 , so that € * is determined by 

/."> /O /O 

— 4- y A - 1 

We now write a', &', o' and /// for a, 6, c and /x, where yd a'&V = /x, ahc, 
since the volumes of the two shells are to be equal. Lastly, putting u + e 3 = w, 
x/a = x'/d 9 &c, we find after some easy reductions that the partial potential of the 
heterogeneous homoeoid at an external point is 

27T/JL L (t) f °° abc chv ■**-./ ax' by' cz 



TT/ Z7T/JL L (t) [ abc cum ,, . / 

y g _ — M ^ 



«-lL (20 J e , Q x ^ \a 3 + w* W + w ' c 3 + w/ ' 

/o /-I /O 1 



uJi = ivR' = w < 



/o /-■) /.-> 

a 3 -f W & 3 + 76' C 3 -j- ?/.'J 



-p. f^ + w ^ 2 . ^ 2 + w rZ 3 c 3 + %o (P 

=- T ; ^ + l^ + lj 2 s + ' ' * + L (n) L (£ + w) 2 3 « + ' * ' 

24. As a first example, fc£ the law of force be the inverse square. We then have 
k = 2 and therefore £ = 0. The partial potential Y ± is then zero and the complete 
potentials, inside and outside, are given by the expressions for V 2 and V 2 ' just written 
down. These agree with those given by Mr. Dyson in the * Quarterly Journal/ 
vol. 25, 1891. 

As a second example, let us find the potential of a thin homoeoid when the law of 
force is the inverse fourth power and the density at (/;, rj } £) is m (£/a 3 ). We then find 
by Arts. (14) and (22), 

\r — ^ 1T / jL ^ t x \ Z -\T _ 27TJJ, f * abc die w 

Vl = T "e m Uy ' Vs = ^~Jo ~qT * 2 Tn. • 

The potential at an internal point is the sum of these two. 
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25. To find the potential of a thin heterogeneous homoeoid when the force varies 
as the inverse K th power of the distance, and k = 2, zero, or any even negative number. 

As the investigation resembles that of the last proposition, except as regards the 
lemma, it is sufficient merely to sketch the argument. 

To avoid confusion as to the sign of k, let us, in the first instance, put k = — k x 
We then have, by Art. 2, 

k y + 1 J J ^ - r 2 

writing as before (Art. 12), 

S 2 d> 
p = $ + r$<j> + r % — + . • • 

V = - --^ 1 j]{<£S Kl+3 + S<£S Kl+4 + . . . }D^ M 

= — *-- 1 f f ^'""^ S„D*da», 
^ + 1 J J L (w — a: x — o) 

where 2 implies summation from k x + 3 to go . We now substitute, as before, 

a r i jt 7 V* v ^ ~y """ jj n^-sz-i ^/ 

to " - IT-'rl - Z L(/)L(7i-27- 1) V Q> 

where j9 = — 2D 2 F, q — ED 3 . 

The lemma by which the integration of each term is effected is as follows. By 
Art. 10 we have 



JJ( a / 2 + /3m 2 + fyn*y +t+1 ~~ ' J \dvj {a 



Following the same notation as before, we let 



A 2 yu, a j/ 2 1 * 1 



+ i> $ + v y + v J Q 



A = A d/c?X + B d/efyi + C d/<fo>. 



We then find, as in Art. 19, 



L(2 S - #) 2tt (- 1)< Aa f , , /<A< ., 1 



where <jr is positive and less than 2s. 

To apply the lemma we put g= n — ^ — 3, 2s — g= n — 2f — l,s+ t — n—f—l. 
These give t = ^(k x + 2) = -| (2 — k) and this requires that the index k of the force 
should be less than 2. 

After completing and summing the same binomial series as before, we have 
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k x + 1 v ' L (n — 2t — 1) J o \clvj 



B»-'-i 



L(»-*-l) Q ' 



where, as in Art. 20, 



X3 3 -2 



By using the theorem of Art. 21, we can perform the operations indicated by A. 
Noticing that the index of the operator (viz , n — 2t — 1) is here less than that of the 
subject (viz., n — t — 1), we find 



V=- 



9. 



irfjb 



>00 



rfV 



/_ lV v l dv ( ~ ) M yr <& ( X 



a 4- ^ 



y 



P + v ' I ' 



M 



K' 






E +1 T 1 



L (0 ' L (t + 1) 2 2 



02 I * * * 



J^t+n ^p» 



L (< + ») L (») 2 3 " 



"J" . . . 



We now write X = a#, &a, and make the necessary changes in the meanings of 
R and T. We thus have 



T 



R 



V 



] f d? dr d 2 ' 

( a + v ) x; 3 + (0 + v ) ^3 + (r + *>) 



^.x 3 



% 



^ 



am* 3 ftvif yvz* 

a + v ft + v 7 + v 



26. We may conveniently express the result in terms of the axes of a, b, c of the 
ellipsoid instead of the squared reciprocals. We find that, when the index k of the 
force is an even integer < 2 or negative, the potential of a thin heterogeneous 
homoeoid, the density at (£ 77, Q being \Jj (£/a 9 r)/b 9 tjc) at an internal point, is 



\T _- ^' 7r / jL 



ft 



1 



foo 




u f du 



[y? 4 



d Y ccbc u* 



\ du] Q x 



Mxjjl 



ax 



b ll 



cz 



\a? -f u 5 3 -f w 6* 3 -f %/' 



where 



1# 
M = — - + 






„D 1 -J- J- J 1 '*" (mD >" ! _l_ 

' 2^ t~ • * • "^ L(« + ») L(») 2 2 » "f" * * •' 



and R, D, and Q x have the meanings given to them in Art. 22. Also t = \ (2 — k). 
At an external point the potential is 



V = 



27T/t 
K-1 



CO 



£\2 



(•it — e 3 )~ f du \ (u — e 2 ) 






abcu 1 



Qi O - e 3 )' 



Mi/> 



&«# 



a 3 -f % 



, &c. ), 



where (x, y, %) are now the coordinates of the external attracted point, so that € 3 is 
defined by 

a 3 + e 3 + ¥ + e 3 + c 3 + e 3 — 
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27. To find the potential of a solid homogeneous ellipsoid at an internal point ivhen 
the latv of force is the inverse K th power and k is even and > 4. 
Referring to Art. 4, we see that when k > 4, 



v = 2(^TTo^T) \\ {uri + ur1 dw - 



Proceeding as in Art. 8, we have 



« i" + V , L (» -/_- 1). /2F\»~y /_ 1 V 
n ^ L (/) L> - 2/) \ E / \KD 



o. / J 



where 25 implies summation fromjf^ to \n. Writing n = k — 4, and remembering 
the integral given in Art. 8, we have 



where 



(a: - 1) L (* - 3) \E/ L (/) 2 a 

= cr — + 6" 77 + <r -—- , 

dor ay dz^ 

and 2 implies summation from y= to \ (k — ■ 4). Thus the series has but two 
terms when the law of force is the inverse sixth, three terms when the inverse eighth, 
and so on. 

28. To find the potential of a homogeneous solid ellipsoid at an internal point when 
the law of force is the inverse k //i and k is even < 4, and may be negative. 

Here 

V = —. =T7 Tn f [(>\ i ~ K + fa*"") do >> 

2 (* — 1) (k — 4) J J v l ' z ' 

r jL±liL — * - L <>-/-1) / _ 2FD« V'-v (ED 3 )/ 

Substituting and using the integral form given in the lemma, Art. 10, we find, after 
summing the same binomial series, 



V 



*•(- iy r _, . id V k <+i 






(* - 1) L (* + 1) J \<fo / Q 
where 

T> — . vou ° 2 vJ3y* uyz* 

a. -{- v /3 -\- v 7 + v 

and 

t = 1 (2 - k). 
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29. Let us next turn our attention to the attraction of a heterogeneous solid 
ellipsoid at an internal point. The density at any point x, y, z being <j> (x, y, z), the 
density at a neighbouring point x + £, y + ??> z + £ is 

Let us describe a sphere of small radius r, whose centre is the point x, y, z. The x 
attraction at the centre of the matter filling this sphere, supposing the law of 
attraction to be the inverse K th and the density to be m£ 9 is seen by an easy 

47T F T* ~ K ~\ r 

integration to be — m I —- except when k = 4. If k is greater than 4, this 

attraction is infinite. If k = 4 the integral is — flog rl and is again infinite. The 

attraction at any internal point of a heterogeneous ellipsoid is therefore in general 
infinite whenever the law of force is such that the index k is greater than 3. We 
shall therefore confine our attention to those laws of force which make k less than 4. 
When the ellipsoid is homogeneous the x attraction of the matter filling the 
elementary sphere is zero at the centre, and there is no restriction on the value of k. 

30. To find the potential at an internal point P of a heterogeneous ellipsoid whose 
density is <j> (x, y, z), ivhen the law of force is the inverse K iJ \ and k is even and less 
than 4. 

Referring to Art. 4 we see that the potential at P of the conical element whose axis 

is in the positive direction of the radius vector is — — — ~^rr , the limits being 

fC "■""* JL J 

r = to r. As in Art. 12, we write 



r 



fy& /ytM 



P=* + r8* + 8»*+...+^8^+.. 



A- K + 



The potential of the positive conical element is therefore — —, t -^ ~^K . 

/c — 1 4 -— /c + m L (m) 

To find the potential of the negative cone, we let R be the distance from P of one of 
its elements Q', taken positively. The density at Q' is found by writing — R for r 
in the above series for p. The potential of the negative cone is, therefore, 



>4 — k + m 



£ — /_ iV^P an( j h ere we mU st write — r 9 for R. 

k — 1 4 — /e + m x ' L (m) ^ 

Taking the two conical elements together, and writing n for 4 — k + w, the 
potential of the whole ellipsoid at P becomes 

- 2 (*.-i)JJ a ^ ^ iT^r+^ri)- 

The integral is to be taken all round the attracted point. The summation repre- 
sented by £ extends from n, = 4 — k to infinity. 
mdcocxcv. — a. 6 c 
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The investigation now proceeds on the same general lines as in Art. 25. By 
Art. 8 we have 

S1±Jjl _ v L(ft-/-1) / _ 2FD 3 V^ 2 /(ED 2 y 

n *L/.L(tt-2/) V } { }} 

where tlie sum extends from jT = to \n or \ (n — 1). Substituting and effecting 
an integration by the use of the lemma of Art. 25, we find, after putting t~\ (2— -k), 



ij 



?y< + r 8 » 8 TC+K - 4 , _ ( - l) lt+t 27r / A^y- 24 - 2 r _, , m /^V -E/(E - E)» JL_ 

n~ L(n + K - 4) L(n - 2t - 2) \ 2 J ) V \ dv) ' L (/) . L (sj * Q ' 



where 



■n 2 o 2 

XV, zn JlL "~r* ~T" '" ' ' — 4" 

A = A — + B-f-' + CI 



> 



<^X ^yC6 ^Z/ 

The summation implied by the 25 does not make a complete binomial series. But 
since the operator contains differential coefficients of the order n — 2t ~ 2 while 
the subject is of the order 25 where 5 = n — /— £ — 1, the terms given by any value 
of f greater than \n are zero. The missing terms of the binomial series are therefore 
zero. Summing the series we find 

v = s <^f (V. *, (f Y p-M^U, i • 

/c -— 1 Jo \av/ L(% -- 2£ -- 2). L(w — # — 1) Q 

where S here implies summation from n = 2£ + 2 to oo . 

To effect this summation we take the series for A^R A given in Art. 21, where g < A. 
Putting <7 = w — 2t — 2, & = w — £_— 1, we find 

L? . LA "" L (2 + 1) L(y) + L (* + 2) L (</ - 2) L (1) 2 3 + * ' ' 

Using Taylob's theorem in the manner described in Art. 21, the result of performing 
this operation on <j> is 

1V+ 1 , Itf +2 T 1 , W + * T 3 1 , 1 , / X . \ 



L(* + l) ' L(* + 2) 1 ' 2 3 » L (* + 3) L (2) '2* ' J r \ * + v 

We now make the changes in the operators described in the article just referred 
to, and the expression for the potential becomes 



(-l)'T _x-7 /'^V 1 tv/tj.7 vx V V 



^=^f>-'*i.R M H 



#2 



\a + -y /3 + v 7 + v, 
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where 

L(t + 1) + L(* + 2) T 2 3 + • ' ' • + L (t + n) L {%) 2 2 * + ' ' ' ' 

avx 2 Bwif <yvz* 



R= 1 - 



T — 

v 3 



1 



' 5 



d 2 cP d* ] 

\(« + v )d*+(f* + v ) & + & + <>)&] 



It is often more convenient to use the axes a, b, c of the ellipsoid instead of the 
squared reciprocals a, /3, y. Let us also suppose that the density at the point x 9 y, z 9 
is \jj (as/a, y/b, z/c). Putting v = 1/w, we have 

*R — 1 _ ^ __ 2/ 2 __ '* a 

a 3 + u & 3 + % c 3 + % 

D _ T^ _ a? + ^ ^ & 3 + ^ J^ c 3 + w d 3 



n»' 



% a 2 dx 2 h 9j dy ?j $ 

Q^ = (a 3 + u) (b 2 + u) (c 2 + w). 
The potential at an internal point of the solid heterogeneous ellipsoid is then 



V 



irabc 



00 



d V v$ ■«*■./ &# % c# 



*du w s r —Mi// 



U - 1 au\ U M : " — lVlt/> — — — j ~ j > 



/c — 1 J o \ #m/ Qi \ $ + w fr 3 + w c 3 + ^ 



~~ L(* ■+ 1) + L(* + 2) T 2 3+,,, + L(HHl)LW 2 2 « + * ' " 

and £ = Jr (2 — k) where k is the index of the law of force. 

As an example, let the law be the inverse square, we then have k = 2 and t = 0. 
The expression then becomes 

+ L(2) 1 2 3 + " " # + L(n + 1) L(n) 2*» + ' * * ' 

which agrees with that given by Mr. Dyson in the ( Quarterly Journal/ 1891. 

31. In Cayley's memoir on Prepotentials we find a simple and elegant formula 
for the potential of an ellipsoid in the form of a definite integral (see ' Phil. Trans./ 
1875, p. 745). Taking the notation of this paper, let V be the potential of a solid 
ellipsoid when the force varies as the inverse K th power of the distance at an external 
point (x, y } z)+ When the density p at an internal point is given by 

6 c 2 
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11 _ "t _ 1 

a? W c 3 



%\n + t 



3 



where a, b, c are the semi-axes of the bounding surface, it is shown that 

where t = |-(/c — 4) and Q 1? € 3 have the meaning given in Art. 31. It is necessary, 
also, that each of the three quantities f + t 9 1 + n 9 1 + t > + -n should be positive. 

When the ellipsoid is homogeneous w + £ = 0, and the restrictions just mentioned 
require that the index k of the force should lie between — f and 4. When the 
index k is greater than 4, the index n + £ of the density p cannot be less than 

i(* — 6). 

Lejeune-Diriohlet (Liouville's 'Journal/ 1839, p. 165) gives a similar formula 
(also mentioned by Cayley) for the attraction of a solid homogeneous ellipsoid, but 
the index k of the law of force is restricted to lie between 2 and 3. 

In the ' Cambridge and Dublin Mathematical Journal/ vol. 3, 1843, Cayley has 
given the value of a multiple integral which, when the number of variables is reduced 
to three, is equivalent to the potential of a homogeneous ellipsoid. Taking as before 
the notation of this paper, and supposing the force to vary as the inverse K th power of 
the distance, the potential at an external point x, y, z becomes 

y/ _ g!gfe A *-l 1 






ifc-i) 2«-» l (0 r (f + o m \_jt— + --£— + 



I (X + a 3 ) 3 ' (X + J 8 )* ' (X + c 3 ) 3 



where * = .£ (« - 4), Q 8 = (X + a s ) (X + &*) (X -f c 2 ), 

^ . y* , ^ -, 



a 2 + X ' 6 2 + \ ' c* + A, 



It is obvious that £ must be an integer, so that k must be an even integer and 
must be greater than 4. 

As the subject of operation is not an integral rational function of x, y 9 z, this 
formula has the disadvantage that the/algebraical work of using the operator is very 

long. 

An elementary proof of this result will be given in Art. 37, by which it will appear 

that it is really a case of Jellet's theorem. 

Another method of finding the potential of a heterogeneous ellipsoid has been 
given by Dr. Ferrers in the ' Quarterly Journal/ 1875. The method requires that 
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the potential V of an ellipsoid, whose density is p == (1 — a£ 2 — firf — y£ 2 )^ should 
be known for all integer values of n. To deduce the potential V 7 , when the density 
at (£, 77, £) is // = $fyfl* he finds a general expression for p f in a series of differential 
coefficients of p with descending powers of n. The value of V is then expressed in a 
corresponding series of differential coefficients of V. The law of force is supposed to 
be the inverse square, but the method applies for other laws of force, provided the 
fundamental expression for V is known. 

32. When the attracted point is external, we may for all laws of force deduce 
the resolved attractions of a homogeneous ellipsoid by the use of Ivory's theorem, 
the attractions at an internal point having been found in Arts. 27, 28. But when we 
know the potential of a thin homoeoid at an external point, we may also express as 
an integral the potential of a solid ellipsoid in which the density p at any point 
(£ V, £) is given by 

p=f(P/a* + 7i*/b* + ?/<?)- 

When the force varies as the inverse square, this problem has in effect been solved 
by Poisson. The potential at an external point (x\ ?/, z') is in that case 



v' = «»&c£{/ 1 (i)-/ 1 K)}^ 



where yi (m) is the integral of/(m), and 

m3 = JTl + vh, + T^ Ql * = {a% + u) {h * + u) (c " + u) - 

Let us turn our attention to the corresponding integral when the law of force is 
the inverse /c th power, where k is greater than the square. 

The potential of a thin homogeneous homoeoid bounded by ellipsoids, whose major 
axes are a and a (I + /a) at an external point, has been found in a finite series in 
Art. 8. Let us use this to find the potential of a thin shell bounded by ellipsoids, 
whose major axes are ma and (m + dm) a. The density of the shell is f(m 2 ), and 
the integral of its potential between the limits m = to I is the potential of the 
heterogeneous solid ellipsoid. 

Referring then to the formula in Art. 8, we unite dmjm for p; ma, mb> mc, for 
a, b, c ; and m 2 a 2 + X, m 2 b 2 + X, m¥ + ^ for a ' 3 > b' l % c ' 2 » We also put X = mhi to 
simplify our results. 

Since the attracted point lies on the confocal.a', b\ c, we have 

T If *Z 

-~ + t~~ + T~, — = m 2 (1). 
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The potential of the solid ellipsoid is therefore 



V' 



where 



1 ale ( 2 V~ 3 



2tt r 

(*-!)(*- 3)J Q 1 \H 



/\m 2 ) mdm \ 1 + — 



1 IMD 



2 3 ac -4 



"~P" « « • 



Q x « = (a 3 + **) (& 2 + w) (« 3 + t*), 



H 



# 



'2 



y 



'3 



*'3 



(a 3 + «i)» "*" (S 3 + uf "*" (c-2 + «)»' 



K 



aV 2 



9 /O 



6V , 



D 



(a 3 + uf 



a" 



dtf* 



(&* + uf d? (e* + uf d 



& 2 $?/ 



>2 



C< 



dfe* 5 



>K* ( «~ 4 >, 



the quantities represented by H, K, D being derived from E' 5 P', A 7 respectively by 
the changes described above. 

Since by differentiating (1) we have 



— 2m 



dm 
du 



x 



'2 



BKloa* 



y 



'2 



(a 2 + uf ' (W + u) 



z >2 

■ 2 T (c 3 + w)» 



H. 



• • (2), 



we may conveniently change the independent variable from m to u, the limits being 
u = oo to w = e 2 where a 2 + ^, b 2 + e 2 , c 2 + e 2 are the squared semi-axes of the 
ellipsoid passing through the attracted point, and confocal with the external surface 
of the given ellipsoidal body. We thus have 



V' = 



2tt 



(*-l) (*-3) 



00 



-afe f(w?)du I 2 V"" 4 f 1 4- -i — — J~ 



> K J (-4) 



where m 2 has the value given by (I). 

This formula expresses the potential of a solid ellipsoid, whose strata of equal 
density are ellipsoids similar to the bounding surface, the law of force being the 
inverse K lh power of the distance where k is not less than the fourth. The series has 
J (k — 2) terms. 

We shalJ now apply this equation to the special cases in which the force varies as 
the inverse fourth, sixth, eighth and tenth powers of the distance. 

33. Suppose the law of force to be the inverse fourth, then k = 4, and the series 
reduces to its first term. The potential of the ellipsoid at an external point is 
therefore 

2ir f 6 * — ale f (m s ) du 



V 
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If the density is uniform, /(m 2 ) = 1, and we fall back on the potential already 
obtained in Art. 5. 

« 

Let us, however, consider what forms of /(nv 1 ) will enable us to effect the 
integrations. Solving (l), Art. 32, we have 

A% 3 + 2Bu -f C 

<mr = 



(a 2 + u) (& 3 + u) (c 3 + u) 
where 

A == x* + y' a + ^ 

2B = {x f * + y'* + s' 2 ) (a 2 + & 2 + c 2 ) - a V 2 - &y 2 - cV 2 

C = a 2 & 2 c 2 (a/ 2 /c& 2 + y /2 /6 2 + 2 2 /c 2 ). 

If we put f(m 2 ) = M/m, where M is some constant, the potential reduces to the 
form 

Writing ^ = l/i; this takes the simple form 

y = — o^&c I77— r~^ 77-* ^i ^ = to 1/e 2 

3 J (A + 2B^ + O 3 ) 4 7 



(A + 2Bs> + CW 8 )* 
where e 2 is given by 



I 1.9, i _2 I „2 i -2 X * 



a S + e 2 l £2 + 6 2 I c 2 + 6 s 

This integral can be easily evaluated, but as the result is rather long it seems 
unnecessary to reproduce it here. 

If the radius vector of the ellipsoid through any internal point R, whose coordinates 
are x 9 y, z, meet the surface in Q, we know that 

where O is the centre. It follows that when the law of force is the inverse fourth 
power and the density at any point R varies as OQ/OR, the potential at any external 
point can be found in finite terms free from all signs of integration. 
This theorem can be generalized. For if we take the density to be 

p = F (m*)jm 

where F is a rational fraction, we arrive at an integral which can be evaluated by the 
use of partial fractions. The potential of the ellipsoid can therefore be found in this 
case also. 
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There is an analogous case where the force varies as the inverse square of the 
distance, and the density of the ellipsoid follows the law p = F (m^)/m. This is 
discussed by Poisson, in the ' Connaissance des Terns' for 1837. After alluding to 
a statement of Jacobi, that the attraction of certain heterogeneous ellipsoids can 
be expressed in finite terms, he shows how this result follows from the formulae for the 
attraction of an ellipsoid which he had previously given in volume 13 of the 
' Memoir es de V Academic des Sciences.' Poisson remarks that the simple case is not 
the homogeneous ellipsoid, but that in which the density varies inversely as m, 

34, Another interesting case is that in which 

M 
(m 3 ) = 



(1 - m 8 )* 
where, as before, M is any constant. Substituting, we have 



= — abc 



f-Md%( x z f a 9 



2 



where 

Q a » = (a s + u) {W + u) {c z + u), 

and e' is the positive root of 

— . + -1— j i— — i (A) 

a* + fl + fl C" + {*> 

If fjb, //, jjl" are the three roots of equation (A), we may obviously write the 
potential in the form 



V ; = — - abc 



^ du - M 



00 



^ (u — /*,)* (^ — niy (to — fji f ) 



i 3 



where ^ = ^ has been chosen to be the positive root. 
We may generalize this result by taking 

p = M(l — m 3 )« + % 
the potential of the ellipsoid is then 

2<tt 7 p 8 — M^A ^ / # V* + * 

V — — abc A — - 1 5— — — «"—** "™" TT~ 

The first of these two cases has some importance, for when the least axis (say c) is 
put zero, the ellipsoid reduces to a homogeneous elliptic disc. The surface density is 
then cr = Mew. The expression for the potential, the law of force being the inverse 
fourth, is then 

y> - fabtr i* ^ \ 1 - 
3 K "Qi I 



X 2 y% $ 1 -f 



($ + ^ & 3 + w w 



CERTAIN LAWS OF FORCE OTHER THAN THE INVERSE SQUARE. 931 



We may verify this result by comparing it with that given by Cayley in the 
c Proceedings of the Mathematical Society/ 1875, for the potential of a homogeneous 
elliptic disc where the force varies as the inverse square of the distance. His 
result is 

V =. 2ab<r ~tt- i 1 



Q 



a? -\- u W -f u 



u 



If Y K be the potential of a disc when the law of force is the inverse K th power, we 

have 

1 dV- 



V „ — 



«-f 1 zdz 



This theorem is proved by Roberts in the ' Quarterly Journal/ 1881. Effecting 
the differentiation indicated, the two values of the potential are seen to agree, 
though the modes of demonstration are totally different. Dr. Ferrers has given an 
extension of Cayley's theorem on the potential of a disc, which, when differentiated, 
also agrees with the results given above, see the ' Quarterly Journal/ 1875. 

35. Let us next consider the potential at an external point when the force varies 
as the inverse sixth power. 

Here k = 6 and the series in Art 32 is reduced to the two first terms. We have 



V 



_2tt r e8 

5.3^ 



ahcf(m*)du ( 2 \ 3 



Qi 



u 



Eu 



IT?/ 

K + ~ DK 

4: . Ji 



where H, K, D have the meanings given in the article referred to. We notice that, 
\ip f be the perpendicular from the centre on the tangent plane to the confocal at the 
attracted point, 



P 



DK = 2^ 



2u du 



1 +^T.+ 



a 3 + u . & 2 + u ' <? -f u 

Q x » = (a 3 + u) {¥ + u) (c 3 + u). 
Substituting in the expression for V' we have 



±dQ 1 
Q x du 



V = £ , f * { - % + 

5.3 i^ { u d Q x 



2p' ^ p't dQ Y 
u 2 Q x uPQj* du 

'2 



die . f (m 2 ) 



where e 3 is the positive root of 



x 



n 



a 3 + € s 



+ 



y 



'2 



Z> 2 +€ 



MOMMM 



2 



'2 



MDCCCXCV. — A, 



3 ' C 3 + 6 3 

6 D 



= 1, 
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If the ellipsoid is homogeneous, we have f(m % ) = 1, and the potential takes the 
simple form 



V 



4tt _________ 

573 e 4 (cfi +75^+^^!^ 



It appears that for all points on the same confocal the potential varies as the 
square on the tangent plane. 

36. When the density is not uniform, we have by an integration by parts 



T77. 



:7T 



5.3 



aoc 



7(1)/ 



«£ 41 



KD 4" 



a'Vc 



f (m 3 ) du 



since m 2 = 1 when tt = e 2 , and by equation (2) of Art. 32 dm*/du 



H = _ \\ v >\ 



We have also written a', V, c for the semi-axes of the confocal which passes through 
the attracted point. 

When the density /(m 2 ) = Mm where M is any constant, we have/' (m 2 ) = M/2m. 
We then find (see Art. 33) 



V = 573 ahc T 



/a 



'■ ' l - ' ■ ■ " ' '** to. ywo 

4 «'7*/„/ • 



r e 2 



cZw 



00 



# (A# + 2Bw + C) J 



The integration can be effected by writing u = l/t\ 

37. The expression for the potential of a homogeneous ellipsoid, when the law of 
force is the inverse sixth power, having been found freed from all signs of inte- 
gration, it follows by Jellett's theorem that the potential for any even power greater 
than 6 can also be expressed in an integrated form (Art. 5). In fact, if Y K be the 
potential of the ellipsoid when the force varies as the inverse K th power, we have when 



V 



Stt ahc 



(* _ 1) L (* - 3) 



A 



i-l 



P 



'2 



o\x 



€ 4 (a 3 + e 3 )* (W + ef (e% + ej 



Writing X for e 2 and remembering the expression for the perpendicular p on the 
tangent plane to the confocal, this is easily seen to agree with Cayleys formulae 
quoted in Art. 31. 

38. When the force varies as the inverse eighth poiver, k = 8, and the series for the 
potential is reduced to three terms. We thus have (Art. 32) 



2tt , 

7.5 



. e 2 



oo 



_^___* / 2 y r K3 h» DK3 w D3R3 

Q_w \i±i(,J [ 4.4. ^ . o . o 



Referring to the values of H, K already written down (Art. 32), we find 



DK 2 =8 



(a* + ^ 



I ... 



+ 4K | 



a 3 J- i(. 



tmJhmm 

| • # * 
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which may be written shortly in the form 



DK 3 = 



8dK 



3 du 



8 J\.v^ 



d 1 
du .Qj* 



where the differentiations with regard to u are partial, the co-ordinates. a?', y\ % of 
the attracted point not being varied. Similarly we have 



D 3 K 2 = 32 Q x 



# 1 
du 2 Q x 8 



Substituting these in the value of V', we find, after some reductions, 



V = — 



87T 



7.5.3 



abc 



f e2 
.(frn?) du 

00 



d 

du 



'2 



P 



du u^Q^ 



If the ellipsoid is homogeneous, the potential at an external point, when the law of 
force is the inverse eighth, is 

87r , t% *L I P' % \ 

77173 ^ fa\fMZj' 

where u 



V = — 



€ 3 . 



When written at length this become^ 



y, 



47T 



7.5.3 



a&c 



JP 



n 



u?Q l 



6 



^ 



+ 



a 2 -+■ ^ 5 2 + u c? -f % 






4p 



2 



03" 



(a 2 + ?0 ? 



y 



2 



+ 



£ 3 



(6 3 + ^) 3 ' (c 3 + u) 



\S 



When the density is not uniform we have, after twice integrating by parts, 



"XT' „_ 



Sir 



7.5.3 



abc 



/(i)p 






/V2 -/*// 






r(m) 



du 



m 2 . 



where w = e 2 in the integrated parts. 

This result takes a very simple form when the density is given by f('M?) = 
The unintegrated term is then zero, and the potential reduces to the two first terms. 

39. When the law of force is the inverse tenth power of the distance the expression 
becomes more complicated. Putting 



U 



P 



'2 



W 



it? d» //\* t / 9 cl 1 4u\ d tp f \* , / 9 <P 1 

\ du Q l Q l J du\u J \ du 2 Q x 



.Qi du* \ u ) 



t d 1 , 4 

4:21 -4— 

du Q l Q x 



P 



u 



/■\4" 



where the differentiations with regard to u are partial, we find, after integrating by 
parts, that the potential of a homogeneous ellipsoid at an external point is 



V 



07T 



9.7,5 3 
6 D 2 



abc u . 
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When the density is not uniform the potential is 



8tt 
~\T / -j 

V - 9T7X3 a0C 



/(l)U + f™* 



40. When the index of the law of force is less than 4, the formula in Art. 32 
becomes inconvenient because the index of K is negative. To find the potential of a 
solid ellipsoid at an external point, we must, in this case, integrate the expression 
given in Art. 26 for the potential of an elementary shell. The integration, however, 
presents great difficulties except in the case in which the law of force is the inverse 
square. The method of proceeding is then as follows : 

To find the potential of a thin homoeoid whose density is xfy ff z h and axes (ma, mb, mc) 
and ((m + dm) a, &c.) we use the formula of Art. 26. Writing dm/m for /x, ma, mb 9 
mc for a, b, c, m 2 co for u 9 and m 2 X for e 2 , we find that the potential of the shell is 

iraoc rf?r £ ' — l ' — — ' * - — - — - - 



A 



L (n) L (n) 2^ \a 3 + wj \6 3 + coj \c 3 + ©/ Q x ? 



where R, D have the meaning given to them in Art. 22, except that o> replaces u, 
and 2 implies summation from w = onwards. 

To deduce the potential of the heterogeneous ellipsoid we integrate this from 

m 2 = to 1. „ i - „ .. 

We notice that 

— -■. . . /YYi "muzz "4"" — -——— «j_ — — 

"*■ — a s + \ ^ If- + \ ^ e? + \ ' 

let jx be the value of X when m = 1, so that a 3 + ju,, ¥ + //-, c 3 + /x are the squared 
semi-axes of that confocal to the boundary of the solid ellipsoid which passes through 
the attracted point, viz., x, y,.z. Then, as ra 3 passes from to 1 , X passes from oo to /*. 
Since R = 1 — m 3 when &> = X, we have 



oo 






,00 



0) 



( ( m a — 1 + B) f /(6)) da. 

J a 



We shall now integrate this equation from m 3 = to 1, i.e., X = co to /x. The 
left side becomes 



^ut, 



w + i/H^-).-lTTT-^) d,B - 



But R and /(w) are both zero for all infinite values of w, hence 

f ^1 /■(„) ^ w = f 1 dm* f (m 3 - 1 + WA°>) d(a - 
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The potential of a solid ellipsoid whose density is xA^z & at an external point is 
therefore 

R» + V»D* / a?x \f I Ihj \ff I ch \* dco 



= Trade 



„, L (n + 1) L (n) 2 2w \a 3 + wj \W + m) \c 2 + «/ Qx # 



This agrees with Mr. Dyson's result, ' Quarterly Journal/ 1891. 

41. When the force varies as an odd power of the distance, the analysis is different 
from that given above. There are two cases, according as the inverse index k is 
negative or positive and greater than unity. 

The former case is the easier, and ultimately depends on integrals resembling those 
used in finding moments of inertia, but of a higher power. It will be convenient to 
take this case first. 

The latter case is the more interesting as we cannot suppose that the force can 
vary as a direct power of the distance except within a limited space. The integrals 
which present themselves in this case are similar to those used in the ordinary theory 
when the force varies as the inverse square. On account of their complexity, it will 
be useful to begin this part of the discussion w r ith the simpler case, in which the 
force varies as the inverse cube, and then to procee^jtojjhe more general law of 
force. 

Though we have here distinguished the cases according as the index k is negative 
or: positive, this is only to obtain finite results. The expression obtained for the 
potential when k is an odd negative integer, may be applied to the case in which k is 
even or odd, positive or negative, but the expression then takes the form of an 
infinite series. 

42. Before proceeding to this formula, it will be useful to notice two other methods 
of arriving at the same potential, 

Mr. W. D. Niven has given in the * Phil.: Trans/ for 1879, a useful formula 
resembling the one used in Art. 43. Let U be any function of the coordinates 
(f 9 y, £) °f a point on an ellipsoid. Taking the integrations over the whole surface, 
we have 

XJp aS = ktrabc I 1 + =-— + — — + . . . 1 U 

where 

V = a 2 dydi 2 + & 2 d*jdyf + c 3 d*ldl\ 

and (£, 7], £) are to be put zero after the differentiations are performed. There is a 
similar expression when the integration extends throughout the volume of the 
ellipsoid, and the strata of equal density are similar ellipsoids. By giving U different 
forms, he obtains the values of some useful integrals. 

To apply this to find the potential of an ellipsoidal" shell at a point (%, y, z) 
for the law of the inverse square, he puts 



U = {(? - xf + ( v - yf +■(£ - zf] 



i 
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He notices that the differentiations with regard to (f 3 77, "£) can be changed into 
differentiations with regard to (x 9 y 9 z) 9 and that, therefore, £ ij, £ may be put equal 
to zero before instead of after differentiation. 

Dr. Hobson has given another method in the ' Proceedings of the London 
Mathematical Society/ vol. 24, 1893. He finds the expansions in series of the 
potentials of ellipsoidal shells, solid ellipsoids and. discs of variable density, the law 
of force being any given function of the distance. Taking the case of a homoeoid, he 
first expands the surface density p in a series of spherical harmonics, say 

pjp = %Y n (x/a 9 y]b 9 zjc) 

where p is the perpendicular on the tangent plane. The potential due to a surface 
distribution on the ellipsoid represented by any one of these terms, when the law of 
force is — <j> (R), is then given by the equation 



Y n <f> (R) p d 



a- 



2 n L (n) 



L (2n ■•*- 1) 



4ira6c(- 1)" t/0m \ t\ \ 1 + ^T7^ + &c> h M a — , b — , c — U(r) 



2 (2tz + 3) 



d ^ d d 



dx dy d: 



where r s == x % + ^ s + ^ R is the distance of the attracted point from the element 
pdcr, and 

A = a 2 d % \dx l + & 3 d 2 /dy 2 + ° 3 d*/dz 2 . 

This series terminates when <^> (R) is an integral rational function of R, and in 
other cases the potential is given by an infinite series. 

43. To find the potential of a thin homogeneous homoeoid at an internal or 
external point P when the force varies as the inverse K th power > where k is an odd 
negative integer. 

Let x, y, z be the coordinates of the attracted point P, £, rj, £ those of any 
elementary mass /xp do- of the homoeoid. Writing 1 — k = 2t to avoid negative 
quantities, we see that the potential of the ellipsoidal shell is 

V = -«- f U d<r[(£ - xf + ( v ~ yf + (£ - *)«]*. 

Let A, R, C, as before, be o?/dx, ^/<%, dJ/cfe ; then the potential is 



V = zyJI jP dae-**-**- *^ + y* + s*)<. 



We shall now expand the exponential and integrate the several terms over the 



CERTAIN LAWS OF FORCE OTHER THAN THE INVERSE SQUARE. 937 

whole surface of the ellipsoid, as, y, z being constant. The terms of the expansion 
with odd exponents need not be considered, as after integration they are clearly zero. 
The integral for any even power is best obtained from the corresponding integral for 
a sphere by the method of projections. We have 

p dor — 5__ _ ^afc v -^ __ l. 9 

J J r L (2s) L (2s + 1) 

where 5 is a positive integer. 

Let A = a* cP/cfc* + & 3 #/% 2 + c* #/cM 

Let the mass ^wfiabc = M. The potential then takes the simple form 






I I 1 + L% + LT5) + • • • + L(^T) + • • • K^ + »* + *^ 



Here 2£ = 1 — k, so that the series is finite if the index of the force is an odd 
negative integer. The expansion is, however, clearly true whether t is integral or 
not, provided only the series is convergent. 

44. To find the potential of a heterogeneous homceoid when k is an odd negative 
integer. 

Let the law of density be p = ^yf^\ Following the same reasoning as before, we 
have 



V = — 



JL 

2t J 



\pdcre' 



A£-B>>-C£/,™2 



(x* + y % + z*)* gift*. 



This is the same as 



V 



ft\\p d<r 



£f{- £f{- so)''- 1 *"'-* ^+ y % + ^- 



Making the same expansion, we have 



— 2t \ dAf 



Cv \t/ 



dB 



d \* a -l - (A% 3 + B 3 6 2 + CVV . 
«C/ L (2s 4- 1) 



When the law of density is p = <f> (£ t?, g), we have 



V = 



4-n-fiabc , I d 
2t ^XdA' 



A. 

dB'' 



d_ 
dC 



^ L(3) T •••>•/■ , 



where r s = a? + y* + 2 2 3 A, B, C stand for -d/dx, d/dy s d'/dz, and 2£ = 1 — *. The 
series is finite when k is an odd negative integer. 

We may also deduce this expression from the series given in Niven's memoir. 
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Putting 

U = £frfC> *{(£ - xf +( v ~ yf + (l ~ *) 3 }, 

the general term of his expansion is 

tirade / # , TO # . d 2 V' T 



L (2i +" 1) \ d? ^ drf^ d£*/ 

Let S = c£/d!£ and A = d/dx. Using the generalized theorem of Leibnitz for the 
differential coefficient of a product, and putting f = 0, after differentiation we have 

the last term being obtained by changing the subject of differentiation from £ to a?. 
Treating yf and £ A in the same way, we have 

<£ (V) U =(- d/dAY(- d/dBy(- d/dCf <f> (AV + B 3 6 2 + CV) ^ (x 3 + f + z%). 

When we replace <j> (V) by V* and \jj by a 2 th power, this becomes virtually identical 
with the general term of the expansion found above for V. 

45. As an example, let the law of force be the direct cube of the distance, and let 
the density be p = £ We have 

~ T 4tirudbc d f . , A% 2 + . . . , , (A? a? + . . -Y A 
4 dk V L(3). T L(5) 

= TT/.a&ci— Arf + T 4 *4 (A¥ + B 9 6 8 +. C*c 8 ) A?- 4 

[L (o) L ^o) 

Remembering the meaning of A, B, C, this is clearly 

V = ~fe- a 3 6c . x (5r» + 8a 3 + 6 8 + c 2 ). 
The potential for a homogeneous homoeoid for the same law of force is 

V = - "f aba {3r* + 4 (a 2 aj 3 -f &V + <***) + 2 (a 2 + 6* + c 3 ) r»} , . ; . 

a constant term being omitted. 

46. We may use the same method to find the potential of a solid ellipsoid when k 
is any odd negative integer. 
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Proceeding as before, we have 



V = - -^ [dv {(£ - xf + (ri- yf + (£- zfy = - ^dv e-**-**- * (a? + y % + z % )t 



where dv is an element of volume. Either by Dikichlet's theorem, or by projections 
from a sphere, we have 



dv 



(A? + B v + C£)* 4*rdbc (A% 3 + BW + C% 3 > 



L (2s) 



2s + 3 L (2s + 1) 



It immediately follows that for a homogeneous solid ellipsoid the potential is 



V = — 



4mabc 



2t 



r 4%2 _|_ (A 3 « 3 -f ) 2 

' • \ \ + -weTT 1 + ""TT^Tir" + • • • f ( x% + ?/ 3 + ^ 



(3). 5 



L (5) . 7 



where A> B, C stand for djdx, d/dy, d/dz as before. 

For a heterogeneous ellipsoid whose density is p = <jf> ( £ ?/, £), we have 



T7 A.irabc . 

V = ~ "27" * 



dA' 









A% 3 -f 
x + -^T^y + . . .\ (x* + y* + z z )'. 



Both these expressions are finite where t is positive, i.e., the index of the power is 
a negative odd integer* 

47. TAis method may be applied to other bodies besides the ellipsoid, provided only 
the value of the fundamental integral 



I 



du(Af-f Br ; + CC) J 



is known for that body. There are many bodies for which this integral is known, such 
as triangles, quadrilaterals, tetrahedra, and so on (see a paper by the author in the 
' Quarterly Journal/ 1886). As an example, consider the tetrahedron, let (£ l9 t^, £ x ), 
(£ 3 , r/ 2 , £ 2 ), &c, be the coordinates of the four corners, let a = A^ + B^ + C£ 1? 
ft — &-%% + B% + C£ 3 , &c, so that a, /?, y, 8 are proportional to the distances of the 
four corners from the plane of reference. The fundamental integral I is then equal 
to the product of the volume of the tetrahedron by the arithmetic mean of the 
homogeneous products of a, /3, y, 8 of the s ih degree. 

By choosing a proper origin and axes the expressions a, /3, y, 8 may be simplified to 
suit any particular case. Using this expression for the fundamental integral the 
potential of a tetrahedron for any direct odd law of force can be deduced from the 
differential coefficients of r 2t . 

48. When the index k is not very large, we can often write down the potential 

MDCOCXCV. — A. 6 K 
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by using the theory of equimomental points. We have already seen that the 
potential of any body is 

V = - I \dv {(g - xf -f ( v - yf + (C~ z)*Y. 

Expanding the t th power, we see that the highest powers of (£, rj, £) which occur 
are in the term independent of x, y, z. These only give a constant term to the 
potential. Neglecting these terms we see that if two bodies are such that l&rf£ u dv 
is the same for both, when f-\-g-\-h is = or < 2t — 1, these bodies have equal 
potentials at all points. Now 2t — 1 is the index of the force taken positively. 
Hence, if two bodies are equimomental for all powers up to the K th , they are also 
equipotential, except for a constant. 

There are many bodies which may be replaced by equimomental points. For 
example, we can show that, up to cubes inclusive, a solid ellipsoid is equimomental to 
six particles placed at the extremities of any set of conjugate diameters, each equal to 
one-tenth of the mass of the ellipsoid together with a seventh particle placed at the 
centre of gravity equal to two-fifths of the mass of the ellipsoid. Placing the six 
particles at the extremities of the axes, we can write down the potential of the 
ellipsoid when the force varies as the cube of the distance. We clearly have 

- 4V = -&M {& + y" + (z - cfY 2 + . . . + f M (a* + f + z*)\ 
therefore, omitting a constant, 

V = - ifeM {5r* + 2 (a 3 + W + c 3 ) r % + 4 (a¥ + &Y + cV)). 

In the same way a tetrahedron, up to cubes inclusive, is equimomental to eight 
points. We collect nine-fortieths of the volume at the centre of gravity of each face 
and one-fortieth at each corner. The potential of the tetrahedron when the force 
varies as the direct cube of the distance can therefore be immediately written down. 
Let the centre of gravity G of the tetrahedron be the origin, and let the axis of z pass 
through the attracted point P. Then GP = z, and (£ l9 7) v ^), (f 3 , %, Q, &c, are the 
coordinates of the four corners„ The potential is easily seen to be 



_ M 



% l ^ . s~c . rt . ~ cm % 



* + ~ % (p + ^ + 3D - ££ s (e + v" + c 2 % 



49. To find the potential of a thin homceoid at an internal point P when the law of 
attraction is the inverse cube. 

Referring to Art. 3, the potential V for an odd inverse /c th power is 



V = -£- r -* ±^ 8 — J)t dm. 

/C — 1 J J 7^ — -Tg 
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When k = 3, this becomes 



Y, = % 



2D' 



do*. 



r l - T 2 



Taking the standard quadratic (3) of Art. 1, we have 



r i - *» 



2D J 



= JF 2 + ^|* = {{odx + £wiy + yws) s + E (aP + ftm + yri>)}\ 



The potential required is, therefore. 



V =- i 



4/* 



f 

I / W CU . 8 » • < • • • . * \xl» 



where ?■' is the radius vector of the auxiliary ellipsoid 

(axX + fyY + yzZf + E (aX a + 0Y* + yZ») = 1 . . . . (2), 
which may also be written in the form 

(1 - fiy* - yz 2 ) aX 3 + . . . + 2/3yyz YZ + . . . = 1 . . . . (3). 



To find the integral we require the following lemma : — 

Lemma. Let r be the radius vector of an ellipsoid ivhose semi-axes are a b\ c f ; 
r" the radius vector of the polar reciprocal with regard to a concentric sphere of 

radius m. Then 

fit / ,» /» 

ab e n .,» , / \ 



abc [ f 

r da* = ~~ A - r" 2 dco 
m 4 J J 



To prove this, we notice that 




v f dxo = dr' da = 




dx dy dz 

x* -f ?y 2 + £ 2 J 



the integration extending throughout the volume of the ellipsoid. Writing x = a'£ 9 
y = 6^, 2 = e'£, and integrating throughout the volume of a sphere of radius unity, 
this becomes 



f 



aVc! dljdr) d% 



r r * 






'7 ' 

— = a o c 



dr dco 



j j j 



c'* cos 3 6 + sin 3 6 (V* cos 3 + a* sin 3 0) 



by changing to polar coordinates. After integrating with regard to r from to 1, it 
becomes obvious that the double intregral remaining is the value of JjV' 2 do* for an 
ellipsoidal surface whose semi-axes are l/a\ 1/6', 1/c'. Introducing the m to preserve 
the dimensions, the theorem follows at once. 

6 E 2 
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50. It- follows from this lemma that the potential of the homoeoid is 

V3— 2m 4 J J 7 a -.-...... (o;, 

where a\ b\ c are the semi-axes of the auxiliary ellipsoid (3). 

This double integral occurs in the ordinary expression for the potential at the 
centre of a homoeoid when the law of force is the inverse square, and its axes are 
m z /a\ w,°°/b' 9 m 2 /c\ Taking its value from that investigation, we have 

s» do = ^ 4 f " du 

where 

Q'*=-(*' + u)(P+u)(y'+u) .(6), 

and a, ft', y are the squared reciprocals of the semi-axes of the auxiliary ellipsoid. 
We, therefore, have 

V3 = Mo! ( 7 )- 

The axes of the auxiliary ellipsoid are given by its discriminating cubic. Let R be 
the leading letter of the cubic, so that the three values of R are a\ /J', y. As we 
want the product represented above by Q'' 2 , we change the cubic into another having S 
for its leading letter by writing R + n = S. The last term of this second cubic 
is the product required. In this way we find 



| (l—/3y%—yz 2 )a-\-u, afixy, ccyxz 

Q'' 2 = | afixy, (l—yz* — ax*)fi-\-u, fiyi/z 

| ayxz, fiyyz, ( 1 — a# 2 — fty 1 ) y+u 



, . (8). 



Substituting either this determinant or the product (6) in (7), we have an expression 

for V 3 . 

If Mbe the mass of the homoeoid, the potential at an internal point when the law of 

force is the inverse cube is given by 

y _ JL f ch l (9) 

We may put this theorem into a geometrical form by remembering that in 
the equation (5) the double integral occurs in finding the potential of a homoeoid 
when the law of force is the inverse square. We have, therefore, the following 
theorem. 
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Let the potential of a thin homoeoid of mass M, and semi-axes a, b, e, be V 3 , when 
the force varies as the inverse cube. Let the potential of a second thin homoeoid 
of mass M" at an internal point be V 2 ", when the force varies as the inverse square. 
If the bounding surface of the second homoeoid be the polar reciprocal with regard to 
a sphere of radius m of the auxiliary ellipsoid of the first, then 

V, V " m 3 



M M" 2abc 

The potential of the homoeoid at an external point P' is the same as that of a 
thin confocal homoeoid passing through P' of equal volume and density at the 
internal corresponding point P. 

51. To discover the geometrical relation of the auxiliary ellipsoid to the homoeoid 
whose potential at P is required, we examine the form (2) of Art. 49. Let OP pro- 
duced cut the homoeoid in C, and let OA, OB, 00, be a set of conjugate diameters. 
Then, if these lines are taken as a new set of axes, the equations of the homoeoid 
and its auxiliary ellipsoid respectively become 

a n y'Z g'2 

OA 3 " + OB* + OC 3 = 1? 
- p i + (^ + ^ + OC,)E=l. 

The two are thus referred to a common set of conjugate diameters. If we take 
A 7 , B 7 , C 7 , in OA, OB, OC, respectively, so that 

/om2 _ /oBy _ _ op 3 , /joey _ j_ f 00 ] 2 _ /2£\ 2 

the ellipsoid having OA 7 , OB', OC' for conjugate diameters is the auxiliary 
ellipsoid. 

52. To find the potential of a homogeneous thin homoeoid at an internal point P, 
when the law of attraction is any inverse odd power greater than 3. 

Putting u = 1/r, the required potential, by Art. 3, becomes 



~\7 __ Hi 

V — — 

ft — 1 , 



*^ 1 K ™ 3 + -263* " 3 do) 
u x — % 2 E 



Referring to the standard quadratic (3) of Art. 1, and putting p = 2F/E, q = 1/ED 2 , 
we have by the theorem already quoted in Art, 8, 

V + u£ =p + npr-*q +...+11 ^^^0^ ) P n ^¥+ • • - 
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This series stops at the first negative power of p and it is necessary that n should 
be greater than zero. If n = 0, the left hand side must be halved. In our case 
we put n = k — 3, and since k is odd the series has \n + 1, i.e., i (k — 1) terms. 

As in Art. 49, we have 

t h - i h _ — ^ + J>3 y _ E ^ , 

where r' is the radius vector of the same auxiliary ellipsoid as that used for the 
inverse cube. In this way we find 

2 (ic — 1) J J A L(/)L(/c-2/-8)VE/ \ED 3 / ' 

where S implies summation from y= to ^ (k — 8). 

In these integrations I, m, n are the direction cosines of the moving radius vector r 9 
and x, y, % are the coordinates of the fixed attracted point P. Let k = ax, /jl = /3y, 
*> = y£, for the sake of brevity, then X, /x, v are constant during the integration. 

Let 

G = r dco (Ik + m ^ + ^) 2<f 

where 25 = /c — 3. As already explained, the integration extends overall positions 
of the moving radius vector round P. 

Remembering that F = Ik + m / A + n v, we see that each term of the series for V 
may be derived from the preceding by a uniform process. Let 

A = a d % jd\* + j3 cP/dp* + y dydv*. 
After some slight numerical reductions we find 



V = 



A* 



*-l2"-* 



1 



IV"" 8 r 4- f±\~\n 1 -u /i.V" fi ^^ 

E/ + \E/ Atjr «-4 + \E/ L (2) (« - 4) (« -- 5) 

i /AV " 3 ~ /A/G L (* zAzi) i 



+ ... 



where jf varies from to J (k -— 3), so that the series for the potential has \ (k — 1) 
terms, k being the inverse index of the law of force. 

"When G has been expressed as a function of X. /x, v, and the axes of the homoeoid 
whose attraction is required, this formula will give the potential at all internal points. 

53. The value of G has been given by the double integral 

G = r' dco (Ik + mijjl + wvf* 
where 2s = k — 3. We shall now show that it can be reduced to a single integral. 
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Let the moving radius vector r' be called OR. The direction cosines of OR, when 
referred to the axes of the attracting homoeoid, are l, m, n ; let its direction cosines 
referred to the axes of the auxiliary ellipsoid he l\ m\ n'. Let the fixed straight line 
whose direction cosines are proportional to X, /x, v, when referred to the axes of the 
homoeoid be called OL. Let its direction cosines, when referred to the axes of the 
auxiliary ellipsoid, be \\ fi\ v. Then 



G = jV dm (cos ROL)*(X a + /x 2 + v 2 Y 



Proceeding as in the lemma (Art. 49) we have 



/ dm (cos ROLf = dr f dm (cos ROLf, 



*l ¥ 



the integration extending throughout the volume of the auxiliary ellipsoid. If x\ y\ z 
be the coordinates of a point on OR referred to the axes of the auxiliary ellipsoid, 
this is the same as 



f 



'd£ dj/clz' /AV 4- fi'/ + v'z'\ 2s 



r 



Putting x r = a'i;, y = b'v), z = c'£, where a\ b', c' are the semi-axes of the auxiliary 
ellipsoid, and transforming back to polar coordinates, we have 



a'b'c 



* » /•> 7 / 



J 



r^dr'do) (\V£ + p'b'ii + vV£) 



(a'f + J'V + ^ +1 



the integration extending throughout the volume of a sphere of radius unity. Putting 
g = r'V, 7) = rm\ £ = rn\ and integrating from r = to 1, we find 



r dm (cos ROL) 3 * = a'-b'c 






Referring to the lemma proved in Art. 10, and putting t = 0, we find that the 
right-hand side of this equation is 



L(2s) 
L(s) L~(s) 2 



sir 

~2* 



*00 



J. 



v~* dv 



r a ,3 V 3 



S'V a 



p 



+ urrz. + 



p* &<•***& 



a'* H~ v b' 2, 4- v g' % + v 



I 



aV-' 



c 



(a* + *>)* (6' 3 4- *>)* (c' 3 + v)* * 



Returning to the notation that a', /3', y (Art. 49) are the squared reciprocals of 
the axes of the auxiliary ellipsoid, and writing v = L/u 9 we have 



r dm (cos ROL) 



1 00 



2s L(2s) 27r r Mt 

Jo 



L (s) L (s) 2* 



M< 



«*<&* 



Q 



/ 9 
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where 



Q'» = ( a ' + u ) (£' + u ) (y + u) 

■v /O /O /D 



■\ '2 '2 '2 

~^-+ >* + - * 

a 7 + tt /3 7 + ^ j + % " 

The separate values of the axes afty of the auxiliary ellipsoid may be found by 
the solution of a cubic, and if this were sufficient, G is here expressed as a single 
integral. But the product Q' has been found in Art. 50 as a determinant without 
requiring the solution of the cubic. We shall now seek for an expression for M of 
the same kind. 

54. Since a', /2', y are the squared reciprocals of the axes of the auxiliary ellipsoid, 

S = (XX + iiY + v z ? + ( aX * + £Y* + r z ') E = 1 . . . . (L), 

it follows that V + u, /5 f + u, y + u & r e the squared reciprocals of the axes of 

S + u (X* + Y* + Z 3 ) = 1 (II.). 

Since these quadrics are concyclic, they are also co-axial, and therefore X', jx,', v are 
the direction cosines of OL referred to the axes of the second quadric. 

It immediately follows that M is the sum of the squares of the projections of the 
axes on the straight line OL and is, therefore, equal to the sum of the squares of the 
projections of any three conjugate diameters on the same straight line. 

Let OP 1? OP 2 be any two conjugate semi-diameters of the section \x + py + vz= , 
and let OP 3 be the third conjugate for the ellipsoid II. Then since OL is per- 
pendicular to the section containing OP l5 OP 3 , the sum of the squares of the pro- 
jections on OL is 

M = OP 3 * (cos P 3 OL) 2 . 

Let X", /a",V' be the direction cosines of OP 8 referred to the axes of the homoeoid. 
Then by equation (II.) 

(«E + v) V 7 (/3E + u) p! f (?E + u) v" 

_ ___ __ p ^ say ^ 



cos PoOL 



/Li V 

xv 7 -f fjLfjjr 4- w ,f 



3— (X» + ^ s + *»)* ? 

1 



OP, 3 



(XX" + //,/a" -f- vv ff ) (XX" -f- fil**" + vv f + />)• 



Substituting these values of OP 3 and cos P 3 OL, we find, after clearing of fractions, 

N 1 



M = 



(U + -«E) (U + j8E) (U + ?E) «f JST A, 3 + JJL? + Z> 3 
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where 

N = X 3 (u + jSE) (u + yE) + p? (u + *E)(u + /3E) + v* (u + «E) (u + /3E). 

If M were written in partial fractions, it must be the same as 



V 3 Zi/ 3 i/ 3 



% + of U + ft U + 7' ' 

as found above. The factors of the denominator must therefore be 

Q'» = (u '+ a) (u + $) (u + y'). 

This is obviously correct since the denominator is the expansion of the deter- 
minant 



Q' 3 = 



u + otE + X 3 Xfx \v 

X[M U + /3E + jJL 9 ' fXV 

\v \lv u + yE + v<l 



■> 



which agrees with the determinant already found for Q' in Art. 50, when we 
substitute..X.^=.aa?,./jt..== /3y, v = yz. 

Finally, substituting in the value of G, we have 

tT " L(s)L(s) t 2 3 * Jo WV .Q' 

where 2s ,= k .— 3 and N, Q' have the values given just above. In this way G has 
been found by a single definite integral in terms of the axes of the attracting 
homoeoid and the coordinates of the attracted point, 

55. The potentials of elliptic rings and discs follow from those obtained for an 
ellipsoidal shell or solid. The transformation may be effected in either of two ways. 
By making the axis of c equal to zero, the homoeoid becomes a disc and the potential 
at any external or internal point may be deduced from the formulae of Art. 12. In 
this case the * surface density at any point £ 77 of the disc is proportional to 
P (1 ■— £ 2 / a2 """ ^V^ 2 )"" 1 where p is the density of the homoeoid at the point £ rj, £. 
In this transformation the law of force remains unaltered so that if the formulae of 
Art. 12 are used, the index must be an even number. 

Another method is to make the axis of c infinite, so that the homoeoid becomes a 

cylinder. Since the potential of an infinite rod, of line density mds, at a distance r is 

1 r£ (* - 2) . r (i) 1 ' ■, ' , . , o 1 1 

~~ZT\ ' vl( _i\ — ^^ 9 we see that the potential of the homoeoid, when the 

index of the aw of force is /c, becomes proportional to that of an elliptic ring with a law 

MDCCCXCV.— A 6 F 
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of force whose index is k — 1, provided the attracted point lies in the plane of the 
ring. In this way the potential of a thin elliptic ring, whose line density at the 
point £, 7), is proportional to pf(£/a, vj/b), where p is the perpendicular from the 
centre on the tangent, may be deduced from the formulae of Art. 12, when the law of 
force varies as an odd power of the distance. 

As these results may be obtained by simple substitution, we shall for the present 
pass them over and consider a theorem which is derived from neither of these 
methods. 

56. When the force varies as the inverse cube of the distance, it is known that the 
level surfaces of a homogeneous elliptic disc are confocal quadrics. The potential at 
any point may then be found in elliptic coordinates by considering the potential at 
some convenient point of each confocal (see a paper by Mr. Webb, 6 Quarterly 
Journal/ vol. 14, 1877). 

We shall now prove that if the law of force is the inverse K th power of the distance, 
and the surface density at f, t? is cr = (1 — ■ ^ja % — ifjWf {K "~ z \ the level surfaces are 
confocal quadrics of the elliptic disc. 

Let P be any point on a confocal, and let the normal at P meet the plane of the 
disc in G, then PG is the axis of the cone enveloping the disc and having its vertex 
at P. Draw two consecutive chords through G ; let one of these be SGS' and let dd 
be the angle between them. Then. PG makes equal angles with the opposite 
generators PS, PS' of the cone. 

Take two elements of area at R, R/ on opposite sides of G, such that PG bisects the 
angle RPR' as well as the angle SPS', let GR = r, GR' = r. These elements attract 

P equally if the surface densities cr, cr at R ? R are such that —^— = ~"~"~~^ — 

Supposing this to be true, the resultant attraction of the two elements, and therefore 

that of the whole disc, will act along the normal PG* 

Because PG bisects both the angles SPS', RPR', we have by elementary 

trigonometry 

1_ 1 1_ _ 1 

~r ~" / ~ GS GS' 

SR . BS' _ sin SPR . sin BPS^ _ sin STB' . sin R'PS SB'.B'S' 

PB 3 sin S . sin W sin S . sin S' PE /3 

The first of these results shows that drjr* = dr'/r'\ Since PR/GR = PR'/GR'? 
the second shows that 

^L _ ^l 3 _ / SB.BS' \»fr- 8 > 

Produce OR, CR' to cut the conic in T, T'. If D be the semi-diameter parallel to 
the chord SGS', we have by a well-known theorem in conies 
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SB . RS' _ _ OR 3 _ _ f _ ^ 

where £ rj are the coordinates of R. There is a corresponding equation for the 
product BR' . R'S'. The necessary relation between or, cr , is therefore satisfied by 
the given law of density. This proof does not require k to be an integer. 

57. To find the potential of the disc, we place the attracted j)oint P on the axis of 
the disc. If h be its distance from the centre C, the potential at P is 

v "-tf-~iJJ (p + 7? 2 + ^r-ir ....... w, 

writing £ = ax, rj = 63/ thus becomes 

^ /c — ijj { C<^ 2 cos 2 + 6 3 sin 2 0) r 2 + A»}*<*- r> * ' ' l '' 

the limits being r = to 1, # = to 27r. This may also be written in the form 

v __ cib Cf rd0 dr (1 - 7£)* {K - Z) / g v 

where a' a = a 2 + ^ ^ 2 = fr 3 + ^ 3 , so that a, V, h are the axes of the confocal through 
the attracted point. When k = 3 the integral is easily found ; we have 

We see from the third form above that 

I. = 1 , _ m -» m\»«-> v, 

the differentiation being effected on the supposition that a, V are constant. 
When the disc is circular, the first form for V„ becomes 



" ~ « - 1 a*~ 8 J J (V 3 + A?)**- 1 ) 



Putting a 3 — r 2 = (a 2 + & 2 ) sin 2 4>, this immediately reduces to 



Vk = ^ __ £ (tan ^- s ^ 



6 r 2 
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where y is the angle subtended by any radius of the disc at the intersection of the 
confocal with the axis of the disc. The integration can be easily performed by a 
formula of reduction. 

58. Conversely, we may inquire whether this is the only law of density which 
makes the level surfaces of the disc to be confocals. Taking the law of force to be 
the inverse K th power of the distance where k> 3, let P be any point on the confocal 
immediately adjacent to the disc. The potential at P is due to the matter of the 
disc in its immediate neighbourhood. Let or be the surface density at this point, 
z the infinitely small distance of P from the disc. The potential at P of an infinite 
plane of this surface density is 

2ttt clr 2it(T 1 

or = 



(k — 1)i*- 1 {k - 1) (k - 3) z K ~ d 

Hence the potential cannot be constant at all points of this confocal unless cr varies 
as 2 K ~ 3 , i.e., or varies as (1 — £ 2 /a 2 — rf/b*)*^"®. 



